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I. INTRODUCTION 


GENERAL INFORMATION 


The present program represents an extension and refinement 
of the previous effort with specific application to the design 
requirements of advanced, small, high-temperature-rise combust- 
ors for aircraft engines in the 2- to 5~pound-per-second (0.91- 
to 2.27-kilogram-per-second) flow range. This program was per- 
formed for the Applied Technology Laboratory, U.S. Army Research 
and Technology Laboratories (AVRADCOM), Ft. Eustis, Virginia by 
the AiResearch Manufacturing Company of Arizona during the period 
July, 1975, to October, 1978. The program is documented in this 


three-volume final report. 


OBJECTIVE 


The primary objective of this program was to further develop 
and validate existing analytical combustor design procedures that 
can be used to significantly shorten the design and development 
cycle of small advanced gas turbine engine combustors. Descrip- 
tions of the combustor analytical models, element tests and model 
validations are presented in Volume I. 


The basic approach of the program consisted of a concen- 
trated analytical treatment of key combustion phenomena affecting 
combustor performance complemented by rig tests. The rig test 
culminated in a complete series of performance mapping to vali- 
date the empirical/analytical combustor design procedure in an 


environment matching an actual operating engine. 


The program was initially comprised of four technical tasks: 


Task I - Analytical-Model Refinement 


4 
! 


ee] 


‘ 


Task Li - Full-Scale Combustor Design, 
Preliminary Tests 


Task III ~ Combustor~Pertormance Mapping 
Task IV - Limited Modification and Retest 


The Task I technical effort is described 
complete description of 


presented in Volume II. 


in Volume I. A 
the Task II and Task III activities is 
The computer codes for combustor design 


that evolved from that effort are fully documented in Volume III 


of this report. ‘The combustor performance goals were achieved in 


Tasks II and III; thus Task IV was cancelled. 


The computer models are based upon the numerical solution of 
the governing aero/thermo equations applicable to 


turbo-= 
propulsion combustor environment, and are, 


therefore, applicable 
for analyzing internal flow field of can, can-annular and annulus 


combustor geometries. Both the inline and reverse-flow combustor 
configurations can be analyzed. 


The cost-eftfectiveness of the empirical/analytical design 
procedure was to be demonstrated by undertaking the design and 
development testing of two full-scale annular combustors based on 


the following enygine/combustor configurations, 


parameters and 
goals: 


l. Engine/Component Configurations, 


@ Annular-combustor contigurations 
e Centrifugal compressor (last stage) 
e First-stage axial turbine 


e Nonregenerative cycle 


Fabrication, and 


if 


Parameters and Goals. 


Engine airflow, W.37 2.87 pounds per second (1.30 kg/s) 
Combustor inlet pressure (P3) « 10 atmospheres 
Compressor efficiency 2 78.4 percent (total-to-static) 
Combustor inlet temperature = 660°F (622°K) 


Combustion efficiency = 99.5 percent (100 percent 
power) = 98.0 percent (5 percent power) 


Combustor pressure loss Po37 Pra = 3 percent 
P 
T3 
Combustor discharge temperature (T sayg) = 2300°F 
(1533°K) 
Maximum pattern factor (PF) <0.23 


T - T 
where PF = A ax fava 
4davg 3 


Average radial temperature profile compatible with 
typical turbine blade requirements 


Maximum radial pattern factor (RPF) .0.075 


where PF = ™4_avg rad max gee avg 

4 avg ~ T; 
"4 avg rad max — peak value of the circumferentially 
averaged radial temperature profile 


e Good light-off/relight capability to 20,000 feet (6091 
meters) altitude and ambient-temperature conditions 
per MIL-E-5007D Paragraph 3.2.5.1 (dated 15 October 

‘ 1973) 


es No visible carbon formation with hot fuel or at high- 
altitude conditions 


e Multifuel capability, including JP-4 and JP-5 


e Fuel contamination tolerance per MIL-E-8593A, Table xX 
with filtration to 10 microns 


e The combined CO and HC exhaust emissions will be suffi- 
ciently low to meet the previously noted combustion 
efficiency goals at 100- and S-percent rated power. 
The NO, LTO emissions level will be at or below the 1979 
EPA NO, standards. The maximum smoke number will be 
below the threshold of the exhaust plume visibility 


e Acceptable component temperature levels and gradients 
to ensure long combustion system life and reliability 


e Reasonable cost and weight 


SUMMARY 


: A complete description of the following six combustor ana- 
lytical models, associated computer codes, and users manuals are 
given in this report. 


e Annulus-flow model 


® Combustor-performance model j 


® Liner-cooling model 


e Transition-liner mixing model 
e Emissions model 
e Fuel-insertion model 


The annulus-flow model is used to compute airflow distribu- 


tion around the combustor liner and pressure drop. The informa- 


tion provided by this model on jet velocities and efflux angles 


is used for specifying the boundary conditions of the internal- 
flow computer models. 


A 3-D reacting recirculating-flow model is used for comput- 
ing internal profiles of velocity components, chemical species, 
and temperature of a given combustor design. Effects of detail- 
design changes can be analytically predicted in regard to combus- 


tion etficiency, exhaust temperature quality, and lean blowout. 


Two-dimensional parabolic programs are used for predicting 
liner-wall-temperature levels, mixing rate in the combustor tran- 
sition liner of 


reverse-flow annuiar combustors, 
emissions. 


and gaseous 
A fuel-insertion model is used to compute mean drop- 
let size and size distribution of pressure atomizers, 
simplex, duplex, 
nozzles. 


including 
and air-assist pressure atomizers, and airblast 
The model is also used to compute spray trajectory and 


evaporation rate of a given nozzle design in a specified flow 
field. 


The use of these models is illustrated in Section IV by a 
worked-out illustration of a simple combustor. 
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The DESCR TET TON 2OF  ANALYTICAT. MODEGS 


The six combustor analytical models are described in this 
section. The relevent computer codes are described in Section 
III. The use of these models is tllustrated by an example in Sec- 
tion IV, 


ANNULUS-FLOW MODEL 


An annulusa-flow model is used to compute pressure lossaesn, 
annulus Mach number and associated air velocity, and airflow daic- 
tribution around the combustor liner. 


A one-dimensional analysis of the plenum annulus is conduc- 
ted based upon the qeneralized one-dimensional continuous flow- 
analysis approach of Shapiro! The analysis considers the effect 
of area change, wall friction, d@raq introduced by inserted 
obstacler such as fuel nozzler and service struts, heat transfer 
from the liner wall, and injection or extraction of air from the 
annulus. The analysis is valid for constant specific heat and 
molecular weiqht, 


Following the approach of Shapiro for a small control volume 
around a point P located at a cistance "X" from the comprensor 
discharge, as shown in Figure 1, the following three working 
relations are obtained for Mach number M, staqnation pressure Pye 
and static temnerature T. 


\ 

Shaptro, Ascher He, "he Pynamics and Thermoevnamics of Com- 
Pressible Fluid Flow, Vol. I", Chapter &®, “he Ronald Press 
Company, New York (1953), 
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Figure l. 
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Control volume around P 


Control Volumes for Annulus-Flow Model. 
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Geo = -ym? 1975 4 ae ax + GF + 2 (1-y) GW 
Po : To H Ww 


=Y PAM 


2 


dT = M (y-1) da + [i-ym?\ fieyei m?\ 47 
ref a ba) 8) & 


_v (y-l) Mm? /af dx + ar = 2y aw 
: 57 PAM Pe 


rae 2 
(y=) (1+ym?) ay] 


Where; 


A = Flow area 

Dy * Mean hydraulic diameter 
f = Coefficient of skin friction 
F = Drag force by inserted obstacles 
T = Static gas temperature 


14 


(2) 


(3) 


a ce nee i 


coe roe ee 


| 


W = Mass flow rate 
y = Injected mass axial velocity/mainstream velocity 
Y = Ratio of specific heats 


The above set of equations is written for each of the con- 
trol volumes, shown schematically in Figure 1, applicable to a 
reverse-flow combustor geometry. Appropriate expressions are 
used for skin friction coefficient and drag introduced by fuel 
nozzles and other obstacles in the flow path. The remaining 
unknown variable dW, which will be negative for the flow through 
various orifices, is calculated by using the following approach. 


The orifice configurations used in a combustor liner can be 
broadly divided into two basic categories. 


e Configurations such as swirlers, primary pipes, and 
venturl sections which are either difficult to handle 
analytically or their flow rates are less affected by 
approach conditions, 


e Liner orifices, including flush port, plunged holes, 
and scooped ports are affected by approach conditions 
and are amenable to analytical approach for predicting 
flow rates, jet velocities, and efflux angles. 


The first type of ports are handled by specifying discharge 
coefficients, whereas the liner orifices are handled by using a 
modified analytical approach described by Gurevice?. The 
Gurevice approach is based upon a 2-D potential flow solution of 
a problem shown schematically in Figure 2. 


For an infinitely long slot of width b, the following three 


relations are obtained for the three unknowns, namely n, 8 and Crs 


2Gurevich, M.1., "Theory of Jets in an Ideal Fluid", Pergamon 
Press, pp 52-59. 


} (4) 
c= cos8/ (1-5) 


a 4 TsinB _ cosh in(t + goat) 
ie oo n n YT - cose 
A 
(1-5, + cos’ R ies (1-5, ) + cogsR (5) 
- a 
: 0-34) cos (1-5,,) = cosB 
C 2 
(1-3, } + (1-4) cose ( -%,)+ (1-2 Joosk i 
a er 
(1-3, oan (-5,)>( “Th ) cogs Ja-35) 
Lh 2 cose (6) 
CD b (55 - T 
Wa2 
Wal Val a Vag 


Where: 
a7L ® annulus width change/upstream width 
b/L = slot width/upstream annulus wiath 
n= annulus upstream tlow rate/slot flow rate, WailWy 


c = upstream annulus velocity/slot velocity, Va’ 
B= jet efflux angle 


Cp = discharge coefficient 


| 

2 

Figure 2. Liner Port Model Schematic. | 
i) 
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These equations are applied to combustor liner orifices hv 


maintaining area similarity through the following relationships: 
: b/, ™ Au/Agg 

tt 

i} 


where Ay s orifice area (2 i p? for circular hole) 


Aaa = effective annulus area with boundarv-layer 


| 

| 

| blockage effects 
‘7 

| 

| 


For a given application, the annulus upstream conditions and 
the static pressure inside the combustor must be apecifiee, with 
the above equations, an orifice can he sized to pass a specified 

, flow rate, or the flow through a specified orifice can be calcue- 
lated. The procedure for each i8 outlined as follows, 

For given values of annulus and orifice flow rates and 
velocities, c can he calculated and then the efflux angle can be 
found from Equation 4. For the special case where all 
flow passes through the orifice, n = 1 and 


annulus 


cos Bp c/2 


If the orifice flow is a neqliqible portion of the annulus 
flow, n approaches infinity, 


cos pre 


After the value of 2 


is obtained, Fauation 5 can he used to 


calculate b/L and Au/Aaa from Fquation 7; then from Fquation 6, 
Cy can be calculated, 
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For the alternate problem, with the orifice specifiece, the 
above procedure is used in an iterative solution starting with an 
estimated flow rate (value of n). The iteration is continued 
until n converqes to a small difference between iterations, 


Such an approach has given qood correlation with measure Cc, 
data of the circular orifices, as shown typically in Fiqure 3. 
For plunged orifices and the metering orifices of film-coolina 
geometries, the approach qave only qualitative aqreement. How= 
ever, by multiplying the computed Cy values by 1.48 and 1.4, and 
by assuming @ equal to 80- and O-deqrees, respectively, for the 
Plunged orifices and the cooling slot, the approach gave good 
agreement with the data, as shown in Figure 3, 


With the above procedure for computing dW appearing in 
Equations 1, 2, and 3, it is now possible to write a set of equa- 
tions for each of the control volumes around the combustor Liner, 
as shown schematically in Figure 1. These equations are solved 
iteratively to compute isothermal combustor-pressure drop. To 
this can be added pressure drop due to heat addition which gives 
combustor total pressure drop. 


3=<D_COMBUSTOR PERFORMANCE MODEL 


A 3-D elliptic, reacting-flow model is used for calculating 
internal flow field of gas turbine combustors. The model solves 
governing equations for the following variables, as described in 
Paragraph 6 below: 


e Axial, radial, and swirl velocity components 
e Specific enthalpy and temperature 
) Turbulence kinetic energy and dissipation rate 


——_ 
Hunter, S. C., K. M. Johansen, H. C. Mongia, and M. P. Wood, 
“Advanced, Small, High-Temperature-Rise Combustor Program, 
Volume IT: Analytical Mode Derivation and Combustor- 
Element Rig Tests (Phases I and I1)", AD778766 (1974). 
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Figure 3, 


and Measured C 


D’ 


A Typical Comparison Between Prediction 


+ Vere ee 


and composite fuel fraction 


e Unburned fuel, CO, 
e Radiation-flux vectors 
® Sprav combustion 


Paragraphs 7 and 8 below give a brief description of the houndary 
conditions and the numerical scheme used for solvina the set of 
nonlinear coupled partial differential equations. 


l. Equations of Continuity, Momentum, and Enthalpy. 


A. Continuity. 
a eft 
div (pv) = M spray (8) 


b. x-Momentum, 


div (p¥u - Were grad u) = > ab - < se (Ut aggdiv v) 


3) . : i Wore er la : so Woge i) Ss 


rar he ‘ 
4 


spray 


{ Cc. y-Momentum. 


ao iv - »- e238 
div (pVvv Vege grad v) a FOr (U ore div V) 
d 


Me , 1 aw Pye ov 
-2F (sap + Flt + Sepray | 
4 


ten me 


a. g-Momentum. 


div (Ww ~ Bees grad w) = - PGP $ oyd Wig div F) 

| + ae OME SB) OS Se Moce © (F Rg > D1 

! +2 3 (wtf GY + avy) - PW + Hatt Gary (11) 
e. Specific enthalpy. 

| div (pV h - oo grad h) = 8) (12) 
| where 

f S, represents the sum of all the enthalpy source terms 


for radiation and spray evaporation 


Definition of variables are: 


H 

] Vo = Net gas velocity vector 

| u,v,w ® Velocity components along x, radial and 
circumferential directions 

x,¥,6 8 Axial, radial and circumferential coordinates 4 


p * Static pressure 


af ‘ 
i h = Static specific enthalpy 4 
u v w 

j S , § , § = Momentum transfer from apray | 
Spray Spray “8PFay to the gas phase u, v and w ~ 4 


momentum equations 
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et 


e 
spray = spray evaporation/combustion rate per unit volume 
< 


{ 

: div (eve) = 4 [2 (eoue) + 2 (cpve) + > Wwe)) 
: p = |Ox reud yy «(re xa fe wo 
| 


div (u grad $) = * en (rug) + x (rug) + st ce $4)] 


The effective viscosity oe, is given by Wore 2 Mg + Uy 


where Uy and WH, are the molecular and turbulent viscosities 
: of the fluid, respectively. 


2% Turbulence Model, 


The turbulent viscosity Me is calculated hv usina ai two- 
equation turbulence model that solves qoverning ecuations for the 
turbulence kinetic energy (k) and the d@issination rate (€). The 
governing equations for k and € are: 


div (pVk - Tk, eft grad k) = G, - pe 


’ (13) | 
div (pVe - Vege G¥Ad €) = (CLG, = Cope) (14) 
where 
2 2 2 
+ a du ov + ow + be 
Gt ule (Gy + Gl) + Gat PD} oe 
2 2 
ow au du dV aw ov w 
+ (x t eae) + (sr t a) + GE + egy 7g) J 
r = , 
k,eff ~ "ert/"k,ef€ , 
: (16) 4 
Vepeff * "ott/c, eff is 
Ue = Cy P ke /e 
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Che Cc 


Ye off 
€,0 
numbers for k and, respectively, 


i? 
are the effective exchange coefficients and Schmidt 


and Cy are constants. Ne ote? le age? Pig oEE and 


The k-€ turbulence model is moderate in complexity and is 
considered to be superior to other models havina a similar deqree 


ee a i ee 


af complexity, This model has been extensively used by many 
investigators and has proved to he adequate in a wide randae of 
flow conditions. More advanced turbulence models, such as those 
based upon the Reynolds-stress modeling approach, are not vet 
fully developed to warrant their use in recirculating Flow-field 
problems as encountered in qan-turbine combustors. In addition, 
such an approach will appreciably increase the computation 
effort. 


Fo ct a” ng mc) a 


Recommended values’ for the constants appearing in the ahove 


ee + 


equations are 


a 


0.09 


hi 


Cy 1.44 


ty = 1.92 


“eoote | 999 


vt v, -¢ is calculated trom 
‘ v, ett 
° Ce alee 
(c, Ce, 


| re Qee 
where k is the vonKarman constant taken to be equal to 0.42. 


fe a. Chemical Species Fquations. 


A two-ntep kinetic scheme ts used as represented hv 


Kauations 18 and 19, 
eke ab 5 +) (oO, + Ny) * x CO + 4 ou,o + n(F + YIN, (18) 


. x 
x CO 3 (0, + n Nj) * x CO, tn (19) 


Daft 


23 


M 
- 
| 
ft 

i. 
4 
) 


‘williams, F. Aw, “Combustion Theorv" 


Stoichiometric oxygqen-to-fuel and oxyaen-to-CO 


first and the second reactions are qiven hv 
x 
i: wee ath 
2. (12 x + y) 
is = 0.571 


Combining Equations 18 ane 19, one obtains overall stoichiometry 
equation as 


The corresponding stoichiometric oxyqen-to-fuel ratio is 


32 (x + %) 
fe—_—_- 
(l2 x 4 y) 


Governing equations for fuel an@ CO mass fractions are: 


div (pVme, ~ Tey gee FLAG Me,) = Mm “ 


boas Rey (21) 
div (Ving - loo, eff grad Mag) a Rag (22) 


where Rey and Rao are rate of oxidation of fuel and CO in accored- 


ance with the combustion model explained in paragraph a. 


mitt 
is the rate of spray evaporation per 


evap’ 
unit volume comnuted in 
accordance with the sprav combustion moeel Cescrirntion in Para- 


graph 5. 


Equations similar to the above are recuired For 05, COs, and 
However, by using Shvap-Zeldavich approximation’, one needs 
to solve only one more equation for composite fuel fraction, 


H,0. 


of u0x 


» Atdcison-VWeslev 
Publishing Company, Inc. (1965). 
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ratios for the 


a 4 
7 +A Ny) * x co,+% HO aon(x 4 a No (20) 


ive ey Peu0x ~ Nguox Perox? = “evap ey 
g~ % 
MASES’ Ghuod  G_ = @, 
m 
where @¢ = Mey 7 — 


, and , are the values of # for air and fuel streams, respec- 


tively. 


Knowing the amount of fuel burned (FB = be uox - Mp.,)+ the 
concentrations of CO., Oo, and HO are given by 


= Xx FB _ 44 
Meoa: ~ 44 Sy) a “co 
Mog = b Mey + Ly Mog + 0-232 - (0.232 +1) Oe. 5y 
18 y FB 


“HO “= 2(12 + y) 
Mass fraction of No is given by 


"No = Lm (oy ' Mey + Mog + "Oy sf 1,0) 


a. Calculation of Reaction Rates. 


Equations are needed for oxidation rates of fuel and CO, 
i.e@., Rey and Rag: The turbulent reactive flow is an area of 
intensive research, and a number of mcedels have been proposed to 
predict burning rate of fuel in turbulent environments. A simple 
model is used in the present study as explained in the following 


paragraphs. 


The rate of oxidation of fuel is determined by the minimum 
of the following three equations: 
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| 


PUAN yoo 
’ Sat tS eee r 
aT me eee =; ss . 3 32a 


7 i 
1 ene! 

oo! 1.5 aw ( =) D8 
Rea, ch Bae “ox Mei? - v aay) 
o 24b 
Reu, turb “Ry » Mey &/k ( ) 

m 
a OX 
Reuox, turb | CR, e iy e/k 

(24c) 


Here, Equation 24a is the rate of fuel oxidation as controlled by 
chemical kinetics. Generally, a bimolecular Arrhenous expression 
is assumed for this reaction. However, the Task I reacting-flow 
mapping data was best correlated by using Kquation 24a in conjunc- 
tion with Equations 246 and 24c. 


Equation 24b is based upon the eddy-breakup model of 
Spalding? and expresses the rate of fuel oxidation as influenced 
by turbulence intensity and scale, and concentration of unburned 
fuel. This model is applicable to premixed flames. Since com- 
bustion in gas-turbine combustors is neither fully premixed nor 
entirely diffusion controlled, Equation 24c is postulated, 
Similar to Equation 24, which determines the rate of fuel oxida- 
tion as controlled by the availability of the oxygen. ‘The con- 
stant iy appears trom stoichiometry of the chemical Equation 18. 
One could use i instead of i, without any loss of generality, 
because the empirical constant “Ry will simply be different. 


The rate of oxidation of CO js similarly the minimum of the 


following three equations: ; 


iy U 
. 2 
ache (2 = 
Reo, ch ~ ‘a & Moo Moy © (x) (25a) 4 
Roo, turb ~ “ye Mag ¢/K (25b) os 
Mon : 
R, =C ep —— 1 /k (25c) ! 
Cox, turb BS ty vol 
*spalding, D. B., “Mixing and Chemcial Reaction in Steady i 
Contined Turbulent Flames", Thirteenth Symposium (Inter- ees 
national) on Combustion, ‘he Combustion Institute, 1971. j 
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4. Radiation Model. 


an A six-flux radiation model based upon the Schuster-Hamaker 
a approximation®’/ is used in the present program. It should be 
noted that, as pointed out by Siddall other flux model approxi- 
mations such as Milne~Eddington and Schuster-Schwarzschild can be 
represented by the same form of flux equations with constants 
being ditferent. Therefore, the user can modify the flux equa- 
tions with relative ease. 


The difterential equations describing the variations of the 
fluxes along the six directions are: 


f(r) = (@ is) ty ban Bt (26) 
a (I a ( ) ~ hoe bad I (27) 
ax sen! ats re a 6 

LSet ps = (ta) tT ee ee (28) 
r dx ry rt r : 67 

4d (ec IL): fa} s) Tt 4 aed “abe 21 (29) 
rar r~ : “Yen r 7 6 

Lg (y )s = (a t 8) I takgteet (30) 
r dd a4+7 | 0+ . 6 
1a ; eg 

¥ an (Ty) (a + gs) Tye oe é I (31) 


SRamaker, H. C., “Radiation and Heat Conduction in Light- i 
Scattering Material", Philips Research Report, Vol. 2, pp 
59-67, 1947. 


Tpatankar, S. V., and D. B. Spalding, "A Computer Model for ‘ 
Three-Dimensional Flow in Furnaces", Fourteenth Symposium bal 
(International) on Combustion, The Combustion Institute, BO 
1973. 


8siddall, R. G., “Flux Methods For the Analysis of Radiant 
Heat Transter", Paper presented at the Fourth Symposium on 
Flames and Industry, 1972. 
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where Tye Tige and Ip, are the fluxes along the positive direc- 
tions of axial, radial and circumferential directions, respective- 


ly; Ivor Tor and Ip. are the corresponding fluxes along the 
negative directions. 


a = absorption coefficient defined as radiation absorbed 
per unit length oe 


8 = scattering coefficient defined as radiation scattered 
per unit length 


E = black body emissive power # ot’ 


where ow is the Stefan-Boltzman constant 


Tat bs T.- t My $ Th. ~ T+ ¥ Ty. 


With the composite-fluxes RS: rR" and R? defined as: 


x oi 
R= F (4+ TW) 
Fae 
R= 9 (Tn, + Dac) 
zi 
BS (Tp, + Ty) 


one can reduce the six first-order flux equations into the fol- 
lowing three second-order equations: 


x 
ax se ae) oe (R*-E) +B (2R* =, Re oo R?) (32) 
zr x 
ae (eee ey on 4 ee) + : (2R*=R*-R?) (33) 
z 
z 
cae Gia yap) = AR oe) Ban? ~ R= Re) (34) 
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Once RR Ro, and Rr? are known, the net radiation fluxes in the 
axial, radial, and circumferential directions, Oo”, Qf, and Q? 
respectively, are given by: 


x 
Q ae Sees ees (ore 
a ar% (35) 
= ata dx 
r 
dq? St Ste. 
2 ar™ 
atadt i 
q* = Tg, ~ Tol 
2 1 aR@ (37 ) 


The contribution of RX, Rr’, and Rr? to the source terms of 
specific enthalpy, Equation 12, is given by: 


(s,) = 2a [(R“ - BE) + 
radiation (RX = B) + (38) 


(R2 = £)} 


Since information on the variations of a and s (with other 
quantities such as concentrations of CO, HO, CO and soot parti- 
cles) is often scarce and unprecise, they have been assumed to be 
uniform. However, variable values of a and s can be incorporated 
in the program with minor modifications. 
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5. Spray Combustion. 


It is very important to predict aerodynamic interaction 
between evaporating/burning sprays and flow field insofar as com- 
bustion efficiency, pattern factor, stability, liner-wal] tem-~ 
perature levels and gradients, smoke and gaseous emissions forma~ 
tion are concerned. For example, the presence of smaller drop- 
lets influence flame stabilization as they provide the main 
source of heat in the recirculation zone. Larger droplets, how- 
ever, escape the recirculation zone and are mainly responsible 
for the smoke formation. Measured data by McCreath and Chigier® 
showed that droplets with initial sizes less than 50 microns were 
evaporated in the recirculation zone. The smaller droplets were 
influenced greatly by the recirculation zone velocity field, 
whereas up to 70 percent of the bigger droplets in the 100-200 
microns range escaped the recirculation zone. Their trajectory 
was not influenced by the flow-field velocity distribution. The 
flow-field influence on evaporation and trajectory of medium-size 
droplets, between 50 and 100 microns, was moderate. 


Combustion characteristics of liquid droplets burning indi- 
vidually are significantly different from those burning collective-~ 
ly in a spray. For example, Beer and his associates! showed 
that the burning-rate constant of monosized droplet arrays was 
about half that of single droplets. There was also significant 
reduction??++4 in drag coefficient, Che as compared to that of a 


eae 
McCreath, C. G. and N. A. Chigier, "Liquid Spray Burning in 
the Wake of a Stabilizer Disc", Fourteenth Symposium (Inter- 
national) on Combustion, The Combustion Institute (1973). 


10Nuruzzaman, A. S. M., A. B. Hedley, and J. M. Beer, "“Com- 
bustion of Monosized Droplet Streams in Self-Supporting 
Flames", Thirteenth Symposium (International) on Combus~ 
tion, The Combustion Institute (1971). 


llonigier, N. A., et. al, "Dynamics of Droplets in Burning and 
Isothermal Sprays", Combustion and Flame, V23 (1974). 
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nonreactive sphere On the other hand, Natarajan? showed that 
Uy Of a burning droplet should be calculated ier a nonreacting 
Sphere at the mean propertics and initial diameter. 


he quasi-steady droplet combustion theory with spherical 
symmetry predicts the burning rate constant K to be independent 
of the surrounding gas pressure, where K = - on det and dy is the 
liquid-droplet diameter. However, it has been found experimen- 
tally l3that K increases with an increase in pressure. Raghunandan 
and Mukundal4 critically evaluated the quasi-steady approxima- 
cion, variable gas phase properties and incomplete combustion as 
related to predictions of burning-rate constant, the flame-to- 
diameter ratio and the tlame temperature. The liquid-phase 
unsteadiness lasts for about 20-25 percent of the total burning 
time. lt was shown by the authors that a good correlation with 
the burning rate data could be obtained by taking thermal con- 


ductivity and Cy as a function of concentration and temperature, 


A majority of the reported work has been concerned with com- 
bustion of single component hydrocarbon tuels. GLimited work has 
been reported, such as Reference 15, for combustion of multi-~ 
component fuel droplets; but, these approaches became quite com- 
plicated for predicting spray combustion of complex fuels like 
jet aviation tuels. 


teNatarajan, R., “Experimental Drag Coetticients tor Evapora- 
ting and Burning Drops at Flevated Pressures", Combustion 
and Flame, V20 (1973). 

\SRush, J. uU, and H. Krier, "Burning ot Fuel Droplets at Pres- 

sures Greater than Atmospheric", Combustion and Flame V22 

(1974), 


14 aghunandan, B.ON., and H. S. Mukunda, "The Problem of 


Liquid Droplet Combustion - Reexamination", Combustion and 
Flame V30 (1977). 


ISehyu, R. Re, Ce Se Chen, G. 0, Gondie and M. M, Elwakil, 
"Multi-Component Heavy Fuel Drop Histories in a High Temper- 
ature Flow Field", Fuel, VSL (1972). 
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Figure 4 pictorially presents the approach that has been 
used for spray combustion in the present program. ‘The spray cone 
is divided into a number of sections or rays. Kach ray has a 
particular x, y, and z direction associated with it, depending on 
the orientation of the fuel nozzle in the combustor and the spray 
cone angle. 


The initial conditions for each droplet are that they have a 
velocity as specified by the program user and a direction corres- 
ponding to the particular ray in question. The total fuel-flow 
rate is currently divided equally among the rays, although it 
would be easily possible to do otherwise. For each ray and for 
each droplet size group, of which five are assumed, the droplet 
trajectories are calculated from a force balance assuming the 
drag on the droplet is for that of a sphere. Heat transfer to the 
droplet is calculated using the coefficient given in Equation 39. 


— 


es 
(140. 3Pr°?Re“) ( a) (39) 


m<K 


h = 2 


OR 


where k is the thermal conductivity of fuel vapor, Re is Reynolds 
number based on relative velocity, and D is the droplet diam- 
eter, Until the droplet reaches the boiling temperature, no 
evaporation is assumed to occur; however, once reached, the evap- 
oration rate is obtained from the burning rate constant k. 
Where: 

$e _g (0) 

° dt 


mee 


o” Pe cP, ln(1+B) 


i : Tee") 
B= mass transfer No. = 1 In lic + a ( ry 0 
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NODES FOR 3-D COMBUSTOR 
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Schematic of Spray in Combustor 


Figure 4. 
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Ho Heat af combustion 


i * Stoichiometric O5-fuel mass ratio 


Latent heat of vaporization 
Pp. = Liquid fuel density 


A, and (CP,) the thermal conductivity and specific heat inside 
the flame zone, are assumed to be 
| 
i 


CP, = CP¢ 
and are evaluated at the average of the hoiling and flame temner- 
atures. These calculations are performed explicitly, with care | 
taken that the time step is sufficiently small, until at least 99 | 
percent of the fuel has evaporated. 

The above procedure recuires knowledge of the fuel and air 
properties. The particular values used in the current program 
are listed below, 


a. Droplet~Size Distribution. 


Group Vol. % of Spray Size Ratio (D/SMD) 
1 0-20 06 i 
2 20-40 0.9 
3 40-60 1.2 
4 60-80 1.5 
5 80-100 2.1 ' 
i 
ri 
;} a 
4 
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b. 


e. 


Liguid-Fuel Densit Po). 
qu ue ensity ( ev 


. : KG 
pp = LO00[PRey + 0.208 - 0.00072 TI “$ 
PR., ® 1.076 for JP4 


1+ aL (1 = 0.67F) 

PReg * Specific gravity of residue at 60°F 
F = Fraction evaporated 
Te = Fuel temperature 

Specific Heat of Liquid Fuel (CP.), 


Molecular Weight of Fuel Vapor. 


Interpolated from table below, 


F MW (JP4) 


93.26 
1 114.60 
3 126.61 
5 138.16 
7 150.59 
9 173.21 
0 204.76 


Thermal. Conductivity of Fuel Vapor (Aged. 
he * 1.729 (A+BT¢) (4/M-°K-Sec) 


where A and B are from the following list: 


| . eis Ret 


MWeuel i B 
50 -6,362E-3 53, 5E+6 
100 ~6.358E-3 49; 1E-6 
150 ~6,284E-3 . 46. 6E-6 
300 -6.010E-3 42,3E-6 
f. Specific Heat of Fuel Vapor (CP. ad 
Vv 
J e 
a = 4183.3 (0.153 + 0.00081T.) (garwp) : 
g. Latent Heat of Vaporization (Pyapl | 
139 
Lyap * 30676.6 (1092.88 ~ 1.87.) (J/KG) ! 
h, 


Boiling Temperature of Fuel (T,). 


° 
Ta 7 A fn P, + B (°K) 


where PY is vapor pressure in pascals, A and B are from 
the following list: 


& Evaporated A B 
0 41.026 -114,000 
10 30.857 41.574 
30 27.348 96.534 
50 23.997 146.567 


The spray calculation procedure is briefly described in the 
following paragraph. Referring to Figure 4, the analysis is done 
for each of the rays selected and their subrays. Details are 
given for the control volume around a point P, where it is shown 
for two typical rays identified as Ray No. 2 and Ray No. 3. With 
each of these rays there might be five subrays, or less, depend- 
ing upon the location of point P, initial droplet sizes, and the 
properties of the field through which the individual drops have 


traveled, Depending upon the direction of the ray and the 
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finite-difference nodal volume, calculations for evaporation/ 

burning are done for a number of subqrid points. The droplets 
| are allowed to exchandqe mass, momentum, and enerqy with the 
5 surrounding gas phase. The net amount of mass, enerav, and 
momentum received by the node P is the sum total of all droplets 
passing through the control volume of P. 


6. Calculation of Gas Temperature. 


With the specific enthalpy and chemical species known, the 
qas temperature is calculated as follows: The specific enthalpy h 
is the summation of the enthalpies of individual species, i.e., 


ee Se = ee eG er ee 


- 
h # : my hy i { 
" } i 
‘ d 
yeom. [h eo aX in i 
ivi 1,0 in ey hh) aT) ' 
ro) 
ye " 
na ; my (Nig + i Sos + -, Cy NAT) : 
e) 
Thus, giving 
hs My my Th, (T*) + Cay ee (47) 


where mie Nye Nit Ch » T*, and T are species mass fraction, 

apecifiec enthalpy, heat of formation at a reference temperature 
Tye isobaric specific heat, aas temperature of the previcus itere 

ation, and the unknown qas temperature, respectively, 


Therefore, from Fquation 47 one obtains the followina 
expression for the gas temperature T: 


h- Sym, h, (T*) 
ah iw Di + — i ee (4 8 ) 
rma Cad (1) 
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The variation of Cp as a function of temperature is taken 


as a fourth-order polynomtal of temperature as qiven in Reference 
16. 


7. 


Finite-Difference Solution of the Equations. 


The numerical solution of the nonlinear, coupled, partial 
differential equations can be obtained by using finite-difference 
methods. A numerical solution of the hydrodynamic equations can 
be obtained by two methods, The earlier approach employed for 2- 
D flows was the so-called streamline-vorticity method, Here 
pressure is replaced from the momentum equations by differentia~ 
tion. Stream function (WY) and vorticity (w) replace the velocity 
components and the pressure, thus requiring solution of only two 
instead of three variables: namely, u, v, and p. The equations 
were solved by a point-by-point successive-substitution pro- 
cedure. Since wy ane w are linked at the boundaries by way of the 
no-slip condition, the w boundary specification could be done a 
number of ways leading to considerably different false diffusion 
levels as recently evaluated by de Vahl, Davis, and Mallinson??, 
In addition, problems were encountered jin obtaining fully con- 
verging solutions with nonuniform qrid@ spacing. Since 1973, 
AiResearch has used a pressure-velocity (primitive variable) 


solution approach, which has the following three advantages over 
the # - w method: 


@ It permits computation of variable density flows where 
p depends upon pressure and temperature, 


16 cordon, S.,and B. J, McBride, "Computer Proqram for Calcula- 
tion of Complex Chemical Equilibrium Compositions, Rocket 
Performance, Incident and Reflected Shocks and Chapman- 
Jouguet Detonations", NASA SP-273 (1971). 


LT de Vehl, Davis, and G. D. Mallison, "An Evaluation of Upwind 
and Central Difference Approximations by a Study of Recir- 
culating Flow", Computers and Fluids (1976). 


e It allows unsteady flows to be calculated as easily as 
steady ones. 


e It works for 3-D flows as well as 2-D flows, whereas 
the #=-w method cannot be easily extended. 


Many primitive variable solution methods have been put for- 
ward by different researchers!®, They vary enormously in com- 
plexity, ease of use, efficiency, and applicability. The 3=<D 
combustor performance code is based on the well-tried SIMPLE 
(Semi-Implicit Method for Pressure-Linked Equations) algorithm 
of Patankar and Spalding as described in Reference 19. The 
features of the computer model includes 


@ Solution of sufficiently general singie form differen- 
tial equations 


e Provision for use with different physical models 

e Use of pressure and velocities as the main hydrodynamic 
variables 

® Use of the pressure correction technique 

e Provision of two coordinate (plane and axisymmetric) 
systems 

e Use of nonuniformly spaced grids 


18 anon., "Proceedings of the Third AIAA Computational Fluid 


Dynamics Conference", Albuquerque, New Mexico, June 27-28, 
1977. 


19D atankar, S. V., "Numerical Prediction of Three-Nimensional 


Flows", in Studies in Convection: Theory, Measurement and 
Application, Volume 1, Edited by B. &. Launder, Academic 
Press, 1975. 


@ Use of staggered grid with attendant minimum truncation 
errors 


e Derivation of finite-difference equations by integrat- 
ing the differential equations over finite control vol- 
umes and thus ensuring mathematical compatibility 

. between the finite difference and the original differ- 
ential-equation formulations 


® Efficient line-by-line tri-diagonal matrix solution of 
the difference equations 


2 Unconditional convergence for all Reynolds numbers 
e Provision for under-relaxation 


A typical grid node spacing is shown in Fiqure 5. Finite- 
difference equations for a node are obtained by inteqrating the 
differential equations over a control volume enclosing a grid 
node. For evaluating the convection and diffusion fluxes through 
a control volume face, a linear variation (in the direction 
normal to the face) of the flow properties is assumed. For other 
purposes, a step-wise variation with discontinuities at control- 
volume boundaries is assumed. Net rate of flow of % into the con- 
trol volume around a node P (Figure %) by convection and dif- 
fusion in the x-direction is 


~ [TL = fye bye + Tyg t+ (Lo £44) Byyl 6 


where 
Ty = Cege,g Ax/>x 
_ UJ 
Ly = m,"/Ay (49) 
Ay 80.5. Wry Pky BY 
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Figure 5.) Typical Grid Spacing and Control Value 
Around a Point P. 
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1 pv = 8, + S) 
equation for the variable $» becomes 


the one-dimensional transport 


Pg PC Sie oy Tee ARON ES Th: et Se op 


X- X= Xt X+ “X+ p p 


wT + (1m Fy) byt by + IT fy by] by, t 8, 

The above equation was derived based on reasonable assumptions. 
However, the linear-profile assumption becomes unacceptable when 
Ey gly y ia larqe compared with Tyy because the weiqhting factor 
(Ty, - fy, Ly.) then becomes neaative, implying an unrealistic 
physical process through which raising the value of Oyy could 
lower the value of >. Therefore, it is assumedc that if the con- 
vective flow rates (L) are larqe compared to the diffusion coef- 
ficients (T), the diffusion across the control-volume face is 
zero and the value of » convection is equal to the value at the 
node on the upwind side of the face. With this assumption, the 


coefficient Tyy - fy byy is replaced hy TE, - fyy Ly 4 where 


My 


Re T (Tyre Cm fxg) bya fyy yy! (50) 


Here [aye Ags a3] stands for the largest of the three cuantities 
Aye Agr and Ay. 


The final finite-~difference ecuation is reeuced to 


Alb @ A 


p'p x+?xe + A 


x-?x- FT Ayatyy t Ayity. + Agyde, + Ap ity + Sy 

(51) 
The solution of the above eauation is obtained bv line-bv-Lline 
relaxation using an efficient tri-diagonal matrix alaqorithm. Av 
this method, for an x-y plane, a traverse along one @irection, 
say the X-direction, is made with ole values for the v.direction 
nodes. Using this solution as the best estimate, the v-Airection 
is then traversed. The same nrocedure is reneatede for ather v-v 
planes. 
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8. Boundary Conditions. 


The specification of the boundary conditions is done in a 
number of ways depending upon the problem. For the left inlet 
boundaries, velocity, density, and turbulence profiles are either 
experimentally known or estimated. The program can handle anv 
specified profiles. For boundaries of the second kind, where 
gradients and not the values of the variables are specified, the 
program uses one of the following two approaches, In the first 
approach, the boundary value is guessed and continually undated 
s0 as to satisfy the given gradient condition. The second 
approach breaks the link through the boundary to all adjoining 
external control volumes by first arranging for the finite- 
difference coefficient connecting the boundary node to an 
internal node to be zero, and then inserting the correct flux at 
the boundary as a false source of diffusion and/or convection tur 
that internal node. 


At the symmetry plane, the convection and diffusion fluxes 
are zero. Therefore, the convection coefficient Cy_ and the 
exchange coefficient (To¢e) are made zero at the axis of sym= 
metry. For the exit plane, information about some of the vari- 
ables is not available. However, since it is the process occur- 
ring in the calculation domain that decides values of the vari- 
ables which the outgoing fluid will carry, there is no need for 
information at such boundaries. These boundaries are simply 
treated by making the boundary Teff equal to zero. The cyclic 
boundary conditions are used for the circumferential direction. 


The near-wall reqion is qiven a special treatment in the 
program. Since the expression for Teff is accurate for turbulent 
flows only, a means is provided for the inclusion of the correct 
shear stresses and other fluxes at the wall. Therefore, the 
nodes next to the wall are assiqned the following values as per 
an empirical wall law: 


! 
\ 


: a ae —a 
yy > 11.5 Faswalte 3 eae (52) 
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Where 6 is the normal distance of the wall from the first inter- 
ior adjacent node. The kinetic energy of turbulence has small 
diffusion near the wall; hence, Nyall for k is set equal to zero. 
Instead of computing I wall for €, it is calculated for the near- 
wall node by assuming a linear variation of the length scale 
giving the following expression: 


LINER COOLING MODEL 


In order to design a durable combustor with conventional 
materials, the liner-wall temperature levels and qradientsa must 
be controlled. Consequently, it is imperative to have a calcula- 
tion procedure that can be universally used for predicting liner 
wall temperatures. The wall temperature at a point P in a com- 
bustor liner (shown schematically in Figure 6) is determined by 
energy balance on a control volume around P, i.e., 


Cy + Ry “ Co + Ra 


where C and R denote the heat transfer rate by convection and 
radiation, respectively. The subscripts H and C correspond to 
the hot side and cold side of the Liner, respectively. 


44 


1 
) 
3 


*sSTepo;; TeoTZATeuy 7o 


uoTzeoOT Ta4dy pue AOSNquSCD Ae[TNUuUY MOTI-SsSTaAsd FO otTzeueuns yo "9 eanhty 


8 : 
f 13GOW NOISSIWA GNV \ f3cGOW ONIXIW 
13GOW JONVWHOIH3d GOS NOLLISNVEL 


7300W ONI1009 YIN ITI 


45 


A 2-D parabolic program is used to compute the hot-side con- 
vection and radiation heat transfer, the marching @irection being 
x, The following expressions are used for calculating the cold- 
side heat-transfer rate. 


-0.2 
Cy = 0.0268(C, G), Ry, (Ty-T) (54) 
1 4 4 
R. = 96 a SP (55) 
c E 4 Bi, (Le. - 7 w "an 
rn Dp &p 


where Chan! Gant and Tn are annulus air specific heat, mass 
velocity, and temperature, respectively. The length Reynolds 
number is based upon x downstream from the cooling-slot metering 
orifices. &, and fo are the liner-wall and the plenum-wall emis- 
sitivites; D, and Dp are the diameters of the liner and plenum, 
reapectively. o and T, are the Stefan-Boltzman constant and the 
liner-wall temperature, respectively. 


One major advantage of using algebraic expressions for the 
colde-side heat-transfer rates is that the appropriate expressions 
can be used for advanced cooling schemes that increase the heat- 
transfer rate from the cold side. Consequently, the cooling 
schemes, such as multiple impingement, extended-surface qeome- 
tries, and chemically-etched surfaces, can be predicted by usina 
the liner cooling model developed in this program. 


Since a 2-D calculation procedure is used for calculating 
the hot-side heat-transfer rates, the model is strictly applic- 
able to either uncooled liners or the liner walls protected by 
cooling films. The user will need to make approximations in pre- 
dicting wall temperatures downstream from discrete radial jets 
such as the primary and secondary jets. 
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The 2-D parabolic program solves the governing equations for 
: the following variables: 


e Streamwise velocity and swirl velocity 


® Turbulence kinetic energy and dissipation model of 
Jones and Launder 2°, 


e Specific enthalpy 


e Unburned fuel, CO, and total fuel appropriate to the 
two-step kinetic scheme. 


e Composite-radiation flux for the two-flux radiation 
model, 
® Five-droplet trajectories 


The governing equations, as reduced from the set of equa~ 
tions presented in paragraph B for parabolic flows, are trans- 
formed to the following generalized form of transport equations 
for the von Mises coordinate system@-, 


9 do. 2 7. 3 
28) + (a + bw) Fh = be to Hey) + ay (56) 
where au ry hr/ (¥Rn¥) (57) 
b= (x thy ~ rothy)/ (Yao ey) (58) 
C= rhu jeg (59) 
(Yer¥y)” 9°54 ete 
wee (¥en vy) (60) 
20 5ones, W. P., and B. E. Launder, "The Calculation of Low- 


Reynolds Number Phenomena with a Two-Equation Model of Tur- 
bulence". ASME Paper 72-HT-20, 1.971. 


2lpatankar, S. V., and D. B. Spalding, “Heat and Mass Transfer 
in Boundary Layers", Intertext Books, London; 1970. 


47 


Here #. is a qeneralized variable, and a, contains the source/ 
sinks and the other terms in the qoverning eauation that do not 
fit in the convection and diffusion terms presented in Equation 
56. Ww, r, and m" denote streamline, radius, and entrainment 
rate across the boundaries W = ve and We Wine The subscripts I 
and E refer to the inner and outer (external) boundaries of the 
domain of interest. 


The numerical scheme used is a variant of the efficient 
numerics of Patankar and Spaldina as described in Reference 21. 
A brief description of how the model is used for predicting 
liner-wall temperatures is qiven in the following paraqranphs. 


Consider a typical combustor liner and its predicted iso- 
thermal lines for an x-y plane; e.q., in line with primary jets, 
such as shown in Figure 7. The 3-D combustor-performance model. 
predicted a revers -flow region near the dome, as shown by hroken 
lines. This particular combustor has three cooling slots on the 
OD liner wall and four cooling slots on the ID wall. The wali- 
cooling model is used for predictina both inner= and outer-liner 
wall temperature with initial conditions for the three stations 
shown, Cefined by the combustor-performance model, The number of 
x-y planes to be solved cCepends upon a particular combustion 
system, If the predicted combustor Internal flow is hiahiv 
three-dimensional, then one may have to solve as manv planes as 
the number of @=-nodes used in the 3-D computation. However, as 
the flow field approaches a two-dimensional approximation, one 
need not analyze more than a few x-y planes to obtain an accurate 
wall-temperature prediction. 


Since the present wall-coolina model is a two-dimensional 
model, it cannot analyze wall reqions near primary, intermeciate, 
and dilution orifices. Conseaquently, it is advisable to restart 
the model for each panel with initial concitions ar aiven bv the 
combustor-performance model. Care shoule he taken in analyzing 
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the primary panel because of the presence of reverse-flow renion. 
There are two possible approaches for analvyzina this section, 
The first approach, the marching reqion lies between the liner 
walls and u >0O with exchange rate specified for 1 > 0 line. In 
the second approach, one can define a boundary edae for which the 
edge conditions for dependent variables, including radiation 
flux, are defined based unon the combustor performance predic- 
tions. Then the wall-cooling model is run separately for the 
inner and outer primary panels. 


In order to get more accurate wall temmnerature predictions, 
it is imperative to know the precise cooling slot exit condi- 
tions. In addition, one must accurately predict the effect of 
the splash-plate thickness on initial mixing between the cold- 
stream and main-combustion qases, The 3-D elliptic code can be 
used with minor modifications to predict the development of the 
jets exiting from the cooling-slot metering orifices. The effect 
of liner-pressure drop, orifice size and spacing, slot-lin 
length, and height on the slot-exit profiles can be analvtically 
predicted, The 3=D elliptic code can also be used to predict the 
effect of the lip thickness on the initial mixing between hot and 
cold streams, 


TRANSITION-LINER MIXING MODEL 


The overall length of a qas turbine engine that employs a 
centrifugal compressor as the last staqe of compression can be 
minimized by using a reverse-flow combustor. When the engine 
uses an axial turbine as its high-pressure turbine, a transition 
liner is needed between the combustor exit and the stator inlet. 
The transition-liner geometry, as shown in Figure 6, is quite 
complicated in that combustor exit flow, with mainly axial 
velocity component flowing from right to left, is bent through a 
180-degree turn in order to flow into the turbine stator. The 
radii of curvature of both the inner- and outer-Liner walls vary 
as a function of distance along the surface. 
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The flow in a practical transition liner is generally 
stream wise with little separation. Small regions of separated 
flows may exist along the ID transition liner surface. But, to 
insure structural durability of the liner, cooling air is 
generally injected in these separated areas so that the local 
liner-temperature levels and gradients are within allowable 
limits. As burner-exit temperatures increase alone with reduc- 
tion in burner length, cooling-film bands might be used for main- 
taining transition liner-wall temperature characteristics at an 
acceptable level. 


A significant fraction of the mixing of hot streaks with the 
cooler combustor gases takes place within the transition liner. 
With the advent of volume-limited turbopropulsion engines, there 
would be need to predict performance of the transition liner in 
regard to exhaust-temperature quality. In addition, the 
transition-liner-wall temperatures must be determined so as to 
estimate the liner life. A 2=<D transition mixing medel was 
therefore developed for this purpose. The program was adopted 
from the wall-cooling model described in the Liner Cooling Model 
paragraph, It solves governing equations for streamwise velocity 
and swirl velocity, specific enthalpy, turbulence kinetic energy, 
and dissipation. The effect of curvature on radial pressure 
gradient is taken into account. Hdéwever, the wressure elliptic 
effects due to streamline curvature have been neglected. 


GASEOUS EMISSIONS MODEL 


Both the combustor-performance model and the wal)-cooling 
model use a simple kinetic scheme in that the combustion process 
is described by two reaction steps, as given by Equations 18 and 
19. Such a scheme is manageable for complex computer codes like 
the combustor-performance model. In addition, the model predicts 
both unburned fuel and CO which accounts for most of the combus- 
tion inefficiency. Since the enaineer is qenerallv interested in 
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eatimating combustion efficiency of a given combustor desiqn, the 
combustor performance model is adequate for the purpose. How- 
ever, when there is a need for estimating the NO, emissions, one 
should preferably use a detailed kinetic scheme so as to predict 
intermediate species (such as O, N, and OH) that are considered 
important for the NO production. A number of calculation pro- 
cedures exist for estimating the NO, levels; e.g., those 
described in References 22 and 23. The approach taken here is 
explained as follows. 


The wall-coollng model was modified to incorporate the fol- 
lowing l6-step kinetic scheme involving 11 chemical species. 


CyHy + (S4) (0, + 3.76 No)—m(xco + 5) HAO + (Fed) (3.76 Maken) 


CO + OH [== CO, + H (62) 
co + 0, a= CO, +0 (63) 
M + CO + 0 == co, +M (64) 
CO + HAO SCO, + Hy (65) 
O, + H, <= OH + OH (f 4) 
OH + Hy [HAO + A (67) 
0, + HIT OH +0 (68) 
o-8 Ho «0H + H (69) 
O + HO =[50H + On (70) ' 
a H+H + M =H, + M (71) , 
O+0+M 50, +M (72) a 


-) eee anborn, J. W., R.~ S. Reynolds, and H. C. Mongia, "A Quasi- 
: Three-Dimensional Calculation Procedure for Predicting the 
Performance and Gaseous Emissions of Gas Turbine 
Combustors", AIAA Paper 76-642, 1976. 


23uosier, S. A., and R. Roberts, “Low-Power Turbopropulsion Poy 
Su Combustor Exhaust Emissions, Volume 3, Analysis", Technical | 
Report AFAPL-TR-73-36, 1974. 
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O+H + Ms=TOH +™M (73) 
| N + OH [ONO + H (74) 
No + Oj NO +N (75) 
0, +N NO +0 (76) 


A one-step global reaction is assumed for oxidation of fuel to 
CO as described by Equation 56. This reaction step is similar 
to the first reaction of the two-step kinetic scheme used in 
the combustor-performance model and the wall-cooling model. 
The reaction step is slightly different from that nroposed by 
Edelmen“4 where his postulation produces Ho instead of H,0, 
assumed here. A set of four reactions are used to describe 
oxidation of CO. Eight steps are used for reactions involving 
Ho, 05 and their dissociation products. Finally, three reac- 
tion steps are used for the NO production. Although the pro- 
gram uses the l6<step kinetic scheme, a more extensive kinetic 
scheme such as that used by Edelman can be incorporated with 
relative ease. 


The fuel-oxidation reaction is controlled by both chemi- 
cal kinetics and turbulence similar to the scheme used in the 
two-step kinetic scheme described in paragraph B3. The 
remaining 15 reaction steps are controlled by chemical kin- 
etics, although the modified eddy-breakup model could be used 
for these reactions also. 


The numerical scheme used in the emission model is 
Slightly different from that used in the wall-cooling model in 
regard to the way the source term 0, of Equation 56 is caleu- 
lated for the 11 chemical species. For each marching step 


24 paelman, R., J. Boccio, and G. Weilerstein, "The Role of Mix- 
ing and Kinetics in Combustion Generator NO,", Paper pre- 
sented at AICHE Symposium on Control of NOY. Emissions in 
Direct Combustion Power Sources, 1973. 
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\ size OXof the parabolic program, a. for the chemical species 
is computed by using the following 1-D equation 


a 
= sama (77) 


~ 


which is obtained by nealecting the cross-stream convection 
and diffusion terms of Equation 56, Equation 77 is solved for 
each of the species by taking a number of steps for the dis- 
tance AX. Typically, 50 steps are used for each AX. With the 
source terms for the species now estimated, Equation 56 is 
then integrated for each of the species over the distance AX. 
Such a modification results in approximately 70-nercent reduc- 
tion in computation time as compared to the numerical scheme 
used in the wall-cooling model. 


FUEL-INSERTION MODEL 


It may be recalled that a sprav-combustion model is used 
in the combustor-performance code, This spray-combustion 
model includes heating, evaporation, and combustion of the 
spray, as well as the sprav trajectories. The code also 
allows for the exchange of mass, momentum, and eneray between 
the spray and the gas phase. Since a complete solution of the 
3-D combustor-performance model takes a lonq computation time, 
on the order of three hours on the Cvber 174, an inexpensive 
calculation procedure was needed for initial selection of the 
fuel-nozzle characteristics. In addition, such a procecure 
would allow approximate evaluation of different nozzle ceslans 
in a flow Field as computed by the combustor-performance 
model. A fuel-insertion model was therefore developed for 
this purpose. 


are eee 


Fuel-droplet evaporation rate and heat-transter rate of 
the droplet are calculated according to the Priem-Heidmann 
model as described briefly in the following paragraphs, 
Vaporization of the droplet mer lpm/see is given by 


mM = Ag K Pap" (78) 
2r_K 
3 L 1/3 1/2 
Nu = 7 = 2 (1 + 0, 

™ Pgape ( 0.3 S, Re ) (79) 

Fs 

YQ S 1 n wo 

Pap PoP yap (80) 


where Ags K, Nuys Puap! tre Pyap! D, Sy» Ra! ana P, are 
droplet-surface area, burning-rate constant, Nusselt number 
for mass transfer, fuel-vapor pressure, droplet radius, fuel- 
vapor density, diffusivity, Schmect number, Reynolds number, 
and surrounding pressure, respectively. : 


Similarly, heat-transfer rate to the liquid surface Oe 
Btu/see is given by 


ay = AGH (MenTy) 2 (81) 
EN. adda Bol 
Uy 2 EO (1 + 0, Pe Ra ) (82) 
”, 
You oe (83) 
etn} 
2 = Me Ce saohh A, 


where h, Nps k, Pye and Cy vap are heat-transfer coefficient, 
5 ¢ 

Nusselt number for heat transfer, thermal conductivity, 

Prandtl number, and isobaric-heat capacity of fuel vapor, 


respectively. 
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i As in the spray combustion model of the combustor perfor- 

Be 

bj mance code, the spray is divided into five discrete droplet 

| sizes. The physical and chemical properties of the jet fuels 

| are varied as a function of the fraction evaporated, as des- 

4 cribed previously in Paragraph 85. 

a The following expressions for the spray SMD are currently 

| 

i incorporated in the code. These can be easily changed by the 

[: user if desired. 
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where: 


We * Fuel flow 
= Fuel viscosity 
v= Fuel kinematic viscosity 
APs “ Fuel-pressure drop 
o = Fuel-surface tension 
S = Fuel-sheet thickness 
Pa Air density 
W. = Air-assist airflow 
Vv = Air-assist air velocity 
Ve" Fuel velocity 
SP. = Primary fuel-pressure drop 
AP. = Secondary fuel-pressure drop 
Wen * Primary fuel flow 
W = Secondary fuel flow 
Pe = Fuel density 
De * Filming diameter 
Wa. = Air-blast airflow rate 


| 
| 
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III. DESCRIPTION OF COMPUTER CODES 
ANNULUS FLOW MODEL 


The annulus-flow model calculates flow conditions around the 
combustor annulus by solving 1l-D fluid-flow equations and pro- 
vides information regarding annulus axial and tangential veloci- 
ties, heat transfer from the liner wall, flow rates, jet veloci- 
ties, jet angles and discharge coefficients of the various liner 
orifices, as well as the overall liner-pressure drop. Coding 
logic is provided so that the user may analyze can, and axial- 
flow and reverse~flow geometries. In addition, options allow the 
program to calculate the pressure drop for a qiven inlet flow 
rate or inlet flow for a given pressure drop. The program will 
also calculate either the flow through a specified orifice row or 
the orifice diameters required to pass a specified-flow rate. 
Finally, a plot, if desired, can be made giving the flow condi- 
tions around the combustor annulus and through the = liner 
orifices, 


The function of the MAIN program (a computer listing has 
been provided in Appendix B) is to call subroutine COMANN, which 
is the main controlling routine, and to perform file manipula- 
tions in the case of an axial-flow geometry. For this qeometry, 
one item (usually not known) is the flow split between inner and 
outer panels. The user inputs essentially two separate cases, 
one for the OD panel and one for the ID. The program will then 
iterate on the flow split until the inner- and outer-panel pres- 
sure drops are equal. While calculations are heing performed on 
one panel, information about the other is stored on scratch 
files, 


Subroutine COMANN performs the iteration loqic and calls the 
other subroutines as required. Iteration on pressure drop or 
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flow rate is performed until the calculated flow through the 
liner orifices agrees to within 0.05 percent of the inlet flow. 
However, if a solution is not obtained in 20 iterations, the pro- 
gram will stop, as errors in the input or high annulus Mach num- 
bers could make convergence difficult. Of particular importance 
is the variable RELAX, defined at card C0.31. This is a relaxa- 
tion parameter used in the convergence logic and has considerable 
influence on the convergence rate. The simple function provided 
haa worked moderately well for various combustor geometries; how- 
ever, for complex designs or high annulus mach numbers it is 
anticipated that the value of RELAX will need to be reduced to 
obtain a solution. 


The names of the remaining subroutines are descriptive of 
their functions. All data cards are read in subroutine INPUT and 
then are printed out in subroutine PINPUT. INLET calculates the 
inlet conditions to the combustor annulus while LENGTH and FLOW 
perform calculations of annulus flow conditions and orifice flow 
conditions, respectively. FLOW, in turn, calls JET and DCOEF 
which calculate the orifice jet velocity and discharge coeffi- 
cient. Some attention to cards DC.34 to DC.57 in DCOEF is war- 
ranted, As there was only qualitative ayreement between the mea- 
sured and calculated diacharge coefficients, a constant multi- 
plier was applied to the calculated values. Line printer output 
is produced in PROUT while the plots are generated in PICTUR and 
BOXES. 


3~D-COMBUSTOR-PERFORMANCE MODEL. 

The 3-D performance model is a three-dimensional 
recirculating-flow program that is capable of analyzing a variety 
of combustor configurations, including can, can-annular, and 
annular. The deck solves for the three velocity components, U, 
Vv, and W, three species concentrations, including UHC and CO, 


turbulence qualities for the K-€ viscosity model, and three 
radiation fluxes, In addition, the use of primitive variables 
makes modifications to the boundary conditions easy, allowing the 
user to analyze complex inlet geometries. Also provided is a 
subroutine for calculating the trajectories and evaporation rates 
of a fuel-nozzle spray. 


Program MAIN (a computer listing has been provided in 
Appendix C) is divided into two basic sections. Up to card 
MA.167, the routine is concerned with reading the input data and 
converting it to the program's internal units which are Systeme 
International (S.I.)}. The input sequence is covered in paragraph 
B of Section IV so only the units will be discussed. Cards MA.7 
to MA.11 are used to define seven arrays which convert lengths 
associated with dimensions and lengths associated with velocity, 
energy, mass, temperature, pressure, and angles respectively. By 
proper specification in the data statements, the user may employ 
those input units that are most convenient, The output units are 
always §.I. From card MA.168 on, MAIN's function is to call the 
other various routines in their proper sequence. 


Subroutine INITIAL performs some preliminary calculations 
(AL.10 to AL.155), prints the input data (AL.156 to AL.258), and 
defines the initial conditions and some of the boundary condi- 
tions on the various arrays (AL.259 on). In section AL.48 
through AL.78, two arrays, JKIN and IKIN, are defined. They 
merely contain flags which indicate the locations of mass injec- 
tion points. Cards AL.261 to AL.272 contain logic for the 
restart option. If Tape 8 from a previous run is saved and then 
made available for use during a subsequent run, the program will 
read the initial and boundary conditions from it. 


Subroutine ALLMOD contains several entry points which per- 
form miscellaneous calculations pertaining, usually, two the 


60 


boundary nodes where modifications to the standard equation are 
in order. The cyclic nature of the boundary conditions in the @ 
or K direction is evident in FMOD as well as limita to the fuel 
and carbon monoxide mass fractions. VELMOD allows the inlet 
swirl velocity to be increased gradually over a number of itera- 
tions and assures that overall continuity is maintained at the 
exit plane, DENMOD makes alterations to the density at the 
boundaries to maintain the correct mass-flow rate. GAMOD speci- 
fies the wall viscosity values as calculated by the wall func- 
tions. SOMAS is used to initialize an array DIVG which is used 
later in the program, The largest entry point SOMOD contains 
logic for modifying the equation coefficients and source terms 
when cooling slots, walls, and droplet evaporation are present. 
Each variable has its own section and accounts for transfer with 
the walls and mass addition from the evaporating fuel. SOMOD2 
deals only with the 2Z-direction radiation equation and is in a 
section alone as the data storaqe is slightly different for this 
variable. 


Subroutine AUX performa the auxilliary calculations for tem- 
perature, density, viscosity, and source terms. Entry DENS uses 
AU.11 to AU.56 to calculate temperature, Cards AU.52 to AU,56 
limit the values calculated in order to account for disassocia- 
tion and early iteration fluctuations, With known temperature, 
density is then determined from AU.57 to AU,108. VISCO obtains 
effective viscosity from turbulent kinetic energy and diasipation 
and calculates y+ for use by the wall function routine, SOURCE 
contains all calculations for source terms with the exception of 
the aforementioned modifications in SOMOD,. Again each variable 
has ita own section, with coding that is quite straightforward 
and requires no explanation. 


Subroutine AUXRAD performs the same function as AUX except 
that it pertains only to the radiation equations. 


SS SS Se 


SPRAY is used to determine the evaporation rate of the fuel- 
nozzle spray. A large section, from SP.106 to SP.269, deals with 
locating the droplet, determining free-stream conditions, and 
handling the situation where the droplet approaches a boundary. 
Next, various fuel- and free-stream properties are evaiuated (to 
SP.292). The drag forces and time step are then determined and 
used to obtain new velocities and location. If the droplet is 
below the boiling temperature, no evaporation occurs (SP.340 to 
SP.347); but, when the boiling temperature is reached, evapora- 
tion rates are calculated, and the appropriate entries to the 
evaporation array (EVAP) are made. Information concerning momen- 
tum changes due to evaporation are also stored in their respec 
tive arrays and later (SP.382 to SP.425) on a seratch file for 
use when the three momentum equations are solved, 


The coefficients for each variable are generated and the 
solution routine called in subroutine STRIDE, First, equations 
for U, V, and W are handled (ST.117 to ST.632), then the pressure 
perturbation (P') is obtained (SP.633 to ST.714) and used to cor- 
rect the velocities (&P.716 to ST.753) so that mass errors are 
minimized, Then, the remaining variables are solved with the 
radiation equations having their own special section (ST,915 to 
ST.937). 


STRAD is a subroutine used in the radiation model which per- 
forms the same function as STRIDE performed for the other vari- 
ables. 


SOLVE provides a solution to the equations generated in 
STRIDE. A full three-dimensional solution would be time consum- 
iny and would require enormous computer storage, Therefore, an 
approximate solution is obtained by "sweeping" through the field 
several times alternately solving along one direction, while 
holding the values in the other two fixed. The variable ICTDMA 
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(UV) at S0O.36 is used to specify the number of such sweeps. As 
the program converges, and the variables assume their final 
values, the sclution becomes more and more accurate. 


\ 
| LINER COOLING MODEL 
{ 


This program is derived from the 2-D parabolic GENMIX pro- 
gram of S.V. Patankar and D.B. Spalding*!, Modifications have 
included the addition of a two-equation viscosity model, two-step 
reaction scheme, two-flux radiation model, plus subroutines for 
calculating wall temperatures and liquid-fuel-evaporation rates, 


| 
| The basic geometry for which this program has been geared is 
continudus inner and outer walls or continuous inner axis of sym- 
| Metry and outer wall, Other situations may be analyzed provided 
the proper internal modifications are made to the code, 
The MAIN program (a computer listing has been provided in 
| Appendix D) handles several functions, including input, estab« 
| lishment of initial profiles, logic for boundary conditions, call- 
ing the additional routines in sequence, and output, The 
| initial section of MAIN (through MA,369) deals with input and 
i initial conditions. Input beging at MA.44 with the case title 
| 


ps SS 


followed by control indices and grid parameters. More computer 
storage has been provided than is required for the six species 
involved in the two-step reaction-scheme, therefore, the extra 
arrays are zeroed, Various other variables are initialized prior 
| to reading the name list at MA.187. ‘The initial profiles are 
“ read from MA.212 to MA.222 and values are assigned to all the 
i arrays from MA.292 to MA.269. The main marching loop (MA.272) 
Ss begins with the calculation of pressure, temperature, and density 
? (through MA.472). STRIDE(1) called at MA.497 calculates the 


2lpatankar, S. Vv. and D. B. Spalding, 
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physical dimension of y from the transformed cross-stream vari- 
able ww, The forward step size is next determined alonq with 
checks for specified X-locations. The boundary conditions are 
established between MA.528 and MA.743, and in this deck can be 
either an inner wall or axis of symmetry and an outer wall. Sec- 
tions dealing with others are bypassed. The wall temperatures of 
the inner and outer walls, if required, are determined by the two 
call statements MA.747 and MA.748, while STRAD (MA.750) is a sub- 
routine used to calculate the radiation flux. Duct geometry and 
pressure gradient occupy the next section MA.755 to MA.857 and 
provides two methods for pressure gradient calculation which are 
selected by IDPDX. When the value is 01, the program uses a 
guess-correction method, whereas for a value of 02, the program 
immediately corrects the velocity and pressure fields if the duct 
area and flow area differ. Entrainment rates are calculated from 
MA.907 to MA,.934 but are not used in any calculations bv this 
code. DROP, a subroutine called at MA.937, calculates the evap- 
oration rates of the liquid-fuel spray. This is followed by 
STRIDE(2), which performs some preliminary calculations needed 
prior to solving the equations. The remainder of MAIN is devoted 
to printout with the exception of STRIDE(3), called at MA.1195 
which actually solves the finite difference equations. 


Subroutine AUX has two parts; the first (through AU,37) 
calculates the effective viscosity from the two equation turbu- 


lence model, while the second computes the source terms for each 
equation. 


As mentioned above, DROP calculates evaporation rates. Note 
that the input data for the fuel nozzle is read at AUS.15. With 
the location of a particular droplet established (AUS,42 to 
AUS.68), properties of the fuel and free stream are determined 
(AUS.76 to AUS.94). Some preliminaries are performed before 
iterative loop AUS.111 to AUS.151 is entered. Calculations are 
performed until the quessed value of the distance the droplet 
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travels agrees with the calculated value. With this distance and 
the droplet velocity known, the time step and evaporation rate 
can be determined. The proper entries in the evaporation array 
EVAP are then made at AUS.173 and AUS.190. 


STRAD performs all calculations relative to the two-flux 
radiation model. Modifications to the source terms due to the 
: presence of a wall are made at GA.130 to GA.136. The central dif- 
i! ference coefficients are then computed and solved using the 
standard tri-diagonal algorithm. 


WTEMP uses an energy balance on the wall to determine the 
wall temperature. Cold-side convection, cold-side radiation, 
hot-side radiation, and hot-side convection are calculated in 
turn, and a Newton-Raphson iteration procedure is employed to 
solve the resulting heat-flux equation. The cold-side velocity 
and temperature are updated at each marching step, accounting for 
the heat transferred to the annulus. 


STRIDE performs the bulk of the numerical calculations and 
has been documented in literature. 


WF is used to evaluate the Couette-flow-equation solutions 
and to obtain wall shear stress and other transfer data. 


PLOTS is a line printer plot routine. 
TRANSITION LINER MIXING MODEL 


This computer code is derived from the 2-D parabolic GENMIX 
program of S.V. Patankar and D.B. Spalding. The primary modi- 
fications include the ability to have a varying step size across 
the grid, since, for a given number of marching steps, the dis- 
tance traveled along the outer-transition liner is considerably 
greater than along the inner, Other modifications include the ae 
addition of K-e viscosity model. | } 

i 
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A computer code that could handle all qeometrical confiqura- 
tions would be greatly increased in size; therefore, this deck 
has been tailored to the geometry of a reverSe-flow annular- 
combustor transition liner. Other configurations can be analyzed 
if the proper modifications are made to the computer code, 


The MAIN program (a computer listins has been provided in 
Appendix E) handles several functions, including input, estab- 
lishment of initial orofiles, loqic for boundary conditions, call- 
ing additional routines in sequence, and output. The initial 
section of MAIN, through card MA.255, contains input and data 
initialization coding. Note that the x and r values of the 
boundaries are read at cards MA.64 and MA.65, but that the value 
of 2, the actual marching direction, is computed at MA.74 and 
MA.77. ISC JE and IPRNT perform the functions their names suqe 
gest, determining which var!ables are solved for and printed. 
Various constants are initialized in the next few cards, MA.157 
igs of some significance since it is here that the name list is 
read, and finally the various profiles are read in and defined. 
Starting at MA.256, the main marching loop begins with the calcu 
lation of pressure, temperature, and density (ha,293 to MA.315). 
STRIDE(1), called at MA.345, is a subroutine which extracts the 
physical cross-stream dimension from the transformed cross-~ 
stream variable, w. The forward step size is calculated in the 
next section (MA.349 through MA.403), which was necessitated by 
having a smaller step size at the inner boundary than the cuter, me 
plus some checks for specified z-locations. Next, the boundary 
conditions are assigned (MA.406 to MA. 479), and since they are 
always walls in the transition liner, only those appropriate sec- 
tions are entered. The actual duct area is determined in section 
MA.481 through MA.505, plus the area required by the flow. 
Should these two not agree, compensation in the pressure qradient 
for the next marching step will be made. The pressure gradient 
is calculated between MA.507 and MA.598. Two methods are pro- 
vided and are selected by IDPDX. When the value is 01, the pro- 
gram usé3 a guess-correction method (MA.539 through MA.543), 
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whereas for a value ot 02, the program immediately corrects tne 
veliccity and pressure fields (MA.521 through MA.536). The pres- 
sure gradient across the grid due to radius-of-curvature effects 
is calculated at MA.579 and incorporated into the axial-pressure 
gradient at MA.595. Statement MA.645 calls AUXO(0), which cal- 
culates the effective viscosity from the K-€ model. Since there 
is no entrainment for the geometry employed, cards MA.646 through 
MA.670 are bypassed. STRIDE(2), called at MA.673, performs some 
preliminaries necessary prior to the equation solution. The rest 
of MAIN is devoted to outputs of various types with the exception 
being MA.902 where STRIDE(3) is called, solving the equations for 
that inarching step. 


Subroutine AUX performs two functions: it calculates effec- 
tive viscosity (up to AU.37) and the source terms for the equa- 
tions (AU.38 on). Subroutine STRIDE performs the bulk of the 
numerical calculations and has been heavily documented in liter- 
ature. Subroutine WF is used to evaluate the Covette-layer- 
equation solutions near a wall, to extract shear stress and other 
transfer data needed in the solution of the equations, and 
finally PLOTS is a line printer plot routine. 


EMISSIONS MODEL 


The emission model is a 2=D parabolic program 
derived from the GENMIX deck of S.V. Patankar and D. B. 
Spalding. The principal modifications include the addition of 
a l6-step reaction scheme and the avility to handle cooling slots 


and radial injection orifices. 


The MAIN program (a computer listing has been provided in 
Appendix F) is concerned with input, establishment of initial 
conditions, logic for boundary conditions calling other sub- 
routines in sequence, and output. The coding is similar to that 
already described for the liner-cooling model; therefore, only 
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those items unique to the emissions model will be discussed. Two 
additional input items are (1) extra specie~input profiles are 
required (MA.213 to MA.217), and (2) data describing the coolina 

slots and radial-injector orifices is read at MA.223 to MA, 275. 

To accompany the additional input, there is also loqic for the 
special boundaries, conditions associated with the cooling slots 

and radial injections MA,562 to MA.618 and MA.659 to MA.715, 
respectively, When the program reaches the edqe of a slot lip, a 
free boundary is assumed until the flow rate of the slot has been 
entrained, A similar procedure is used for radial-injection 
orifices where the boundary is assumed to be a porous wall. All 
the other subroutines perform the same functions as described in 
the Liner-Cooling Model paragraph; however, an additional sub- 
routine, AUXS, has been added which calculates the specie source 
terms and writes them on a scratch file. AUXS solves the same equa- 


are omitted. The equations are solved many times using a step 
size considerably smaller than the main program. In this manner, 
an estimation of the change in the specic value for the larger 
Main program marching step is obtained from which an averaqe 
source term over the interval can he calculated. This is then 


tions as STRIDE except that cross-stream convection and diffusion | 
I 


used when the complete equations are solved in SYTRIDK(3). The 
first section of AUXS (up to AUS.68) performs data initialization 
followed by the calculation of forward and backward rate con- 
stants (AUS.70 to AUS.95). The main locp is entered next where 
first the source terms and then the derivatives are determined, 


Note that the rate expression for the qlobal fuel reaction con- 
tains the effect of turbulence (AUS.125) and that the kinetic 
source term contains a number of variables each raised ta a 
respective power (AUS.123 and AUS.124). These powers, HFU, ERO, 
etc., are read in through the name list in the MAIN proqram. A 
step size, such that the species values do not chanae exces- 
sively, is then selected using the variables TERM] and TRRM?, 
also part of the nama list. Examination of AUS.223 shows that 
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the maximum change allowed during the marching step is the laraer | 
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of TERM2 and TERMI times the upstream specie value, The Foaua- 
tions are then solved and the process is repeated until a dis- 
tance equal to the main program step size, DX, has been 
traversed. ISMAX, which is also part of the name list in the MAIN 
program, is a limit on the maximum number of these steps. The 
average source term over the interval is then calculated and 
stored on a scratch file for later use, 


FUEL INSERTION MODEL 


The function of the fuel-insertion model is to determine the 
evaporation rates and trajeautories of a fuel spray in a two- 
dimensional flow field. Information concerning the fuel nozzle 
and flow field are read in, and from these the program calculates 
the fuel SMD and the trajectories, Tf desired, a plot of the 
droplet paths can be made and the evaporation rates saved for use 
in other programs, 


INJECT] (a computer listing has been provided in Appendix G) 
is the main program and controls inpnt, output, and the other 
Subroutines, Up to INJ.67 several data statements initialize 
some fuel and air propertien needed in the evaporatian calcula- 
tions. Input data is read to INJ.119 followed by some prelimin- 
ary calculations, Additional input is read at INJ.187 to INJ.195 
if a nonuniform flow field is specified. Calculations of SMD for 
the particular fuel nozzle type selected begin at INJ.206 con- 
tinuing to INJ,310. The next section entered (INJ.311 through 
INJ.435) loops over the five droplet sizes, calculating, in turn, 
their trajectories. The remainder of INJECT! provides output and 
plotting. Subroutine FEVAPC is used to save the evaporation 
rates for later use, if desired, while AIRPRP interpolates the 2- 
D nonuniform flow field to obtain the free-stream conditions, 
Subroutine EVAP performs the majority of the calculations, 
including the force balances on the droplets sq that their 
trajectories and the evaporation rates can be calculated, 
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IV. ILLUSTRATIONS 


ANNULUS FLOW MODEL 


Figure 8 shows the annulus-loss model example geometry. It 
is a simple reverse-flow combustor with inner- and outer-panel 
cooling slots, plunged primary orifices and two dome inlets. For 
the purposes of analysis, the annulus was divided into elements 
denoted by the dashed lines, These divisions correspond to 
places at which mass was extracted or where the annulus was of 
irregular shape. The elements can be length-type for which skin- 
friction losses and heat transfer, etc., are calculated or flow- 
type for which mass extraction is calculated. Thus, the length 
element (2) is the annulus section between the inlet and 
and flow element|3 lis the OD cooling slot. 


Inspection of the input sheets, Figures 9 through 11, shows 
that Card 1 contains the case title and specifies some control 
parameters. The annulus-inlet conditions are specified on Card 2 
along with the total number of elements used and the inlet- 
element number. Card 3 is for internal liner flow only and is 
omitted, whereas Card 4 specifies various constants. It now 
remains to describe the plenum and liner shape and the various 
orifices. This is done on the second input sheet with explan- 
ations of the various items given in Figure 11. Since, for this 
case, the flow spiit between the various orifices was known and 
not the orifice size, the flow-element cards are of the fixed- 
flow ratio type (FF). Had the orifice size been known instead, 
the fixed diameter type (FD) should have been used. If the Cy of 
a particular orifice row is known, it may be specified as has 
been done for flow elements 8 and 10. Even though these are 
actually annular slots in the dome, the program will still cal- 
culate an orifice diameter, which is, of course, meaningless; 
however, the effective and geometric areas are also provided in 
the output from which the correct annular-slot height can easily 
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TITLE - MUST HAVE FOR EACH CASE 78 79 BO 


ANNULUS INLET FLOW CARD 
NEL NELI ; 
1 345678 11 Wh 2, PTL 4, TTL 4, BETAl,)DP/P 6) 7 


2 fos Joc lon? liote. J20 Jor | ——] 


DOME INLET FLOW CARD (REQUIRED FOR INTERNAL LINER FLOW ONLY) 


NEEI oy P rT TA 
12 676 11" 93 PT oan 4 BETA 5) 61 72 


Tie Sa Cee a A Da ee ll Meare 


CONSTANTS CARD 
L 2 PRICS) BLKF 3, TANS 4, CDB 51 SK 6S +, TREAX 


{le es a ee iss lee 


1 Title - Run ident appears on printed output and plot 


*e 1 some (or all) holes tixed, inlet flow fixed, 
iterate to get pressure drop. 


= 2 some (or all) holes fixed, pressure drop fixed 
iterate to get inlet flow (input Wl is first guess). 


j= O no plot 
ZPIC in 1 plot drawn 


= Q output printed after converged solution 
rIpsuaG ie L output printed after each iteration) Use only to de- 
= 2 output printed after each element | bugg 
non-convergence 


NEL 2 total number of alement stations (card 2 = only) 
NELI = element ID no. (NELEM) at annulus or dome inlet 2. and 
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3 
Wy and W = air flow at inlet stations, lb/ sec 
PT, and PT = inlet total presaure, PSIA 


TT, and TT = inlet total temperature, R 
Beta,, Beta # swirl angle at iniat 


DP, P = total presaure drop, PSIA, PSIA (first guess if ITER = 1) 
card 2 only 


4 \" 0, smooth wall friction factor 
FRIC ~l, no wall friction 
= roughness factor for rough walls 
BLKF = annulus effective area factor (= .83 for fully developed 


turbine raat 
TANS = tangent of flow separation spread angle (.1 recommended) 


te ca rt are Ce a em 


Figure 9. Annulus Loss Model Input Sheet (Sheet 1 of 2) | 
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CDB = drag coefficient of struts across annulus (l1~1.2 RECM) 
' SK = ratio of air specific heats, 
i RG = &ir gaa constant 
fis = 0. for reverse flow annular or can combustors 
- ITREAX)= 1, for axial flow annular. First data set is for OD 
: panel. Program expects a second set for ID panel 


e CASE TERMINATION 


After last card of case: 


: © In Column 1, Column 2 blank - case repeated with changes, 
g next card ia title card followed by cards with changes 
i from previous run. 


oo «In Columns 1 and 2, next card is EOF to quit or new title 
ecard followed by all cards for complete new cage. 


se oS 


Figure 9. Annulus Loss Model Input Sheet (Sheet 2 of 2) 
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Figure 10. Sample Work Sheet for Program 117. 
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ALL NUMBERS MUST HAVE 
NELEM DECIMAL, POINTS 


N HOLES 


N HOLES 


Flow panauge ia divided into length (LL) and flow (F) elements, element numbers, NELEM, 
oan be in arbitrary ordar, i.e., LO, L, 3, 4, 16, 30. The cards are atacked in order 
from inlet 4o laet P because numbers are arbitrary, 4a new alement can bo inserted without 
renumbering other cards. 


? flu _ LENGTH EGEMENTS All Dimensions in Inches 


First & card ia annulus inlet (cL = 0) = ror voth external and internal cases. 
QSNQTH ELEMENT 

a tnturnal and external flow cases must be run 

aaparataly 
For internal casee, and card ia dome inlet (LI) and 

LC cards are used with T RISE « 4T due to combustion 

in thhe elument (no & carda) 
For uoxternal cases usa only L Cards 


SLUNUM 
XL» XP = X COORD to end of element L « Liner 
RL, RP » Radius to und of vlement PP » Planum 
(ror internal flow XP, RP = Op, XL, 


RLo« ID) 
CL = Length of element (optional) 
TLIN «= Mean apts temp. over CL, °R 

tf = 0 then TLIN & TTL 
BLKI = . 


COMMUN ORIGIN 
+ tt ALL OLEMENTS 
RANT Area . qt aie 
DAKI = Width of strut 


18 @ 80, STRAIGHT 
VINE 1B CALS. 


ELORDEUCH FLO Sh uMeht 
“arda.ara inserted between L cards at pointe where flow is extractud, Flow cunditions 
into F vlem, are those from upsetroam [. elem, 


Types: FP = Fixad Flow Ratlo, W/W 
FD + Fixed Orf. Diam, (W/Wl fe First Guess) 


FQ» Bleed flow (not included in Liner flow) 
FI - tnternal flow elem. (input co thea. elementa ie obtalnuod from an external 
flow solution) 


WW) © Orifice Plow/Inlet Flow 
NHOLES - No. of Orifices 
NHTYP : Hole Ba hee (For Cp) 


t 
. Plush dole, Thin Wal} Chott. BAAS em 
2. Plunged Hole mene’ 
3. Cooling skirt 4a ee us 
4. Flush Hole, Thick Wall 


aonnae eene nice 
5. Cd Input 
6. Rectangular Hole 

ecb = DISCH Coefficient (NHTYP « § Only) 


a ae 
DHOLES «© tlole Diam. (Por FD Only) 

T SEP + 1, Separation is Reattached 

VIET = Jet Paceien Velocity, FPS 

ANGT «= Jet Injection Angle, Dag. 

ABETA © Swirl Angle in Annulus Outeide FI Elem. 


Figure 11. Program 117 Input Data Sheet 


Input Format for Element Cards 
Sheet 2, 
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be obtained. Other irregular shaped orifices, such as slotted- 
dilution holes, for example, can be handled in a similar manner. 
4 Card deck termination is specified by the double periods. 


The output of the program is shown in Figures 12 through 14. 

‘ Fiqure 12 is merely a printout of the input data. Figure 13 gives 
information of conditions in the annulus at the various stations 
the user has specified. Figure 14 has information regarding the 
orifices. Note that the program has calculated orifice-hole 

| diameters and that the flow split is the same as was specified in 4 
the input. Overall output parameters such as pressure drop, cor- \ 4 

rected flow, etc., are also given. The numbers under the heading 

INPUT FOR SUB BOXES= are associated with the plot subroutine. 


COMBUSTOR PERFORMANCE MODEL 


The combustor geometry for the 3-D combustor-performance 
model example is shown in Figure 15. The reverse-flow annular 
liner has been divided into a qrid network consisting of 30 nodes 
in the axial or X direction, 19 nodes in the radial or y direc- 
j tion, and 13 nodes in the tangential or @ directions. The 
| 


a 


decision of how large a @ segment to analyze is depending upon 

where radial planes of symmetry can be found, since the program 

assumes that there are cyclic boundary conditions in the @ direc- 

tion, Therefore, any mass leaving the system along the K-13 

plane is assumed to reenter the system through the K = 1 plane and 

vice versa. For this example, uniform grid spacing has been 
2 used, although this is obviously not required. 


Cem awe ae a te Se ee ce er 


p | The completed input sheets for the case are shown in Figure 
nn 16. Additional input information can be found on the input sheet 
forms located in Appendix A. Cards 1 and 2 are titles used for 
printout and case identification, The grid size (30 by 19 by 13) 
has been entered on Card 3 along with indicators for axisymmetric 
geometry, K-« viscosity model, kinetic and turbulence controlled 
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Figure 16. 3-D Combustor Performance Model (1 or 6) 
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Figure 16. 3=-D Combustor Performance Model (3 of 6) 
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Figure 16. 3-D Combustor Performance Model (4 of 6) 
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Figure 16. 3-D Combustor Performance Model (5 of 6) 
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Figure 16. 3=-D Combustor Performance Model (6 of 6) 
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combustion, relative pressure, and radiation to be included, 
respectively. The use of relative pressure, rather than absolute 
pressure, merely means that the pressure field is referenced to 
one particular grid node by subtracting that pressure value from 
the entire field. On Card 4, the uSer-selected-input units have 
been specified. The particular units employed are combinations 
of convenient ones and are specified by data statements in the 
main program, at lines MA.9 to MA.11. The user should redefine 
the arrays for his own input units or use S.I. by specifying IU = 
01. Specifying the perfect gas-density calculation, not printing 
the initial values, the inner boundary as a wall and that this is 
not a restart, complete the information on Card 4. In order to 
restart a case, Tape 8 must be saved from the previous run and he 
made available to the program. Card 5 indicates that all the 
variables are to be solved. P' is the pressure perturbation used 
in the solution algorithm, h is stagnation enthalpy, and Fue Fie 
and Fy are the three radiation fluxes. Card 6 shows that all the 
variables will receive one sweep of the solution routine per 
iteration except P', which will have six. The solution of the 
equations in early iterations is not exact, but this does not. 
present a problem since the coefficients are changing from itera- 
tion to iteration anyway. However, it is beneficial to the con- 
vergence rate to obtain a more accurate solution to the P' equa~ 
tion. The variables to be printed are specified on Card 7. The 
06 merely indicates that every 6th X-plane will be printed. The 
relaxation parameters specified on Card 8 have been used suc- 
cessfully on a wide range of problems and probably need not be 
altered, The Prandtl numbers and coordinates, Cards 9 through 
13, are fairly se) Ff explanatory except that the first field, Card 
12, is the radius of the inner boundary and that values of y(2) to 
Y(MPL) are measured from that line, Cards 14 through 16 are used 
to specify an inclined wall at the inlet. Since there is none for 
this case, the values of IWEI and JWIO are set equal to 2 and MP1 
(i.e., 19) respectively, and Cards 15 and 16 are omitted, Were 
there an inclined wall, Fiqure 17 Gives the necessary information 
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| i ACTUAL WALL 
i SMUCAVES~ 
7 ite 
16 
iS IWE! = 08, 
14 JWIO = 09 
je Co masa 
43 SA eS me 5, 6, 5, 6, 8,8, 8, 8 
46 SH ee ce eee 
Pete cette JWLO = 8, 9, 12, 12, 
: 16, 16, 19, 19 
8 tab ip 
? | 
6 
5 
4 


THE ACTUAL WALL IS SIMULATED BY STAIRSTEPS ALONG THE MIDPOINTS BETWEEN 
NODES. JWIO,IS THE JNODE WHERE THE INCLINED WALL STARTS, i.6,,08, AND IWEL IS 
THE I-NODE WHERE THE INCLINED WALL ENDS, i.e,08, IWLI IS THE 1ST I-NODE INSIDE 
THE SIMULATED WALL AT A PARTICULAR J-LOCATION, ELEVEN VALUES ARE RE- 
QUIRED SINCE THIS IS THE NUMBER OF J-NODES INCLUDED IN THE INCLINED WALL, 
J= 9-18, AND ARE MARKED BY. JWLO 1S THE 1ST J-NODE OUTSIDE THE SIMULATED 
WALL AT A PARTICULAR I-LOCATION. EIGHT VALUES ARE REQUIRED SINCE THIS IS 
THE NUMBER OF I-NODES INCLUDED IN THE INCLINED WALL, | = 1-8, AND ARE 
MARKED BY oon 


Figure 17. Information Necessary to Describe an Inclined Wall. 
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required to describe it. Cards 17 through 20 need no further 
explanation other than that given in input sheet forms. A total 
of 150 iteration steps, along with no intermediate printout, have 
been specified in the first two fields of Card 21. The basic pro- 
gram contains provisions for only one dome inlet, specified by 
JSW1 and JSW2; however, the example problem has two. The second 
inlet was handled by internal modifications to the program in the 
Subroutine ALLMOD, a process required on nearly every combustor 
analyzed, as a code which makes fcrovisions for all possible con= 
figurations would be extremely complex and lengthy. Note that 
each primary orifice is simulated by three nodes, making the 
total of vadial-injection points equal to six. Inlet conditions 
for the dome are provided on Card 22, while Cards 23 and 24 
describe the fuel nozzle. A back angle (/3) of zero would have the 
nozzle spraying in a purely tangential and in an increasing #¢ 
direction. Positive B would have the spray cone rotated toward 
| 
| 
| 
| 
| 


the dome; therefore, a value of +90 degrees has the spray cone 
directed axially along the combustor as is required. Zero down 
angle (8) was needed for this example. A positive value would, 
however, rotate the spray cone toward the inner wall. For qeom- 
etries that have the entire spray cone in the calculation domain, 
such as this annular one, THETAL and THETA2 never change from 0 
and 360 degrees or their equivalent in whatever angular input 
units are selected, Cards 25 throuch 30 are used to describe the 
cooling slots,of which two were specified (NUINJ = 2 on Card 21). 
Note that even though the slots surround I-node 16, the I loca- 
tion is specified as one greater (or 17) due to internal conven- 
jences in the program. Similarly, Cards 31 through 38 specify 
the radial injection points. Due to the sign convention on V- 
velocity (positive is in the positive y-direction), injection 
points on the outer wall have a negative V-velocity component, as 
shown on the last three fields of Card 34. 


ee a i coe a 


90 


ek em a ave 


2p ea 


The output of the 3-D model is illustrated in Figure 18 (7 
sheets). Sheet 1 is a printout of some input data, including 
fuel and air flow rates, injection velocities, and other impor- 
tant quantities. Note that the output units here are S.I. 
Sheets 2 and 3 show the u-velocities at @ = 6 deq (in line with 
the primary holes) after 150 iterations. The total error in 
mass, i.e., the sum of the mass error in all the grid cells, was 
1.5 percent of the total flow, which was deemed accurate enough 
for this example. The two dome inlets J = 2-3 and J ® 17-18 can 
be seen at I # 2, The u-evelocities are calculated for slightly 
different control volumes than the other variables. The U- 
velocity printed at I = 8, for example is actually the U-velocity 
that occurs at the cell boundary between I = 7 and I # 8. This 
displacement results in the inlet U-velocity being stored at I = 
2 rather than I= 1. A small recirculation zone exists behind the 
dome and is terminated by the primary orifices at I = 10, 11, and 
12. Small recirculation regions are also evident behind each jet 
at I = 14, while the presence of the two cooling slots can be seen 
at I # 16 and 17. Sheets 3 and 4 of Fiqure 18 show the fuel mass 
fractions for the same @ plane. The figure also shows the 
extremely rich regions just behind the dome at I = 2-4 near the 
fuel nozzle. By the time the exit plane is reached, some 
unburned fuel still remains. Temperature is shown in Sheets 5 
and 6. The fuel-rich region behind the dome (seen in Fiqure 15) 
exists at a relatively low temperature as the amount of oxyqen is 
very limited. Farther down the combustor where the primary jets 
have recirculated (I = 6,7) one sees temperatures closer to 
stoichiometric. The primary jets penetrate to the combustor cen- 
terline as evidenced by the temperature profile at I = 12. These 
jets produce a colder core with hot regions on either side that 
extend well past the cooling slots at I = 16. These slots provide 
a cool film that extends to the exit of the combustor. Sheets 6 
and 7 show the evaporation rate of liquid fuel. One can clearly 
see the two sides of the spray cone and that some fuel is imping- 
ing on the wall at I = 9. Practically all the fuel has evaporated 
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by 
| 
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by To = 18 with a few larger droplets persinting at ot o< 20, where 


they entered trom neighboring K-planes. 


LINER-COOLING MODEL, 


The Liner-cooling model wan used to predict the wall temper- 
atures downstream of the cooling aloets in the combustor shown to 
Figure 15, Starting with the profiles at the exit of the slot, an 
predicted by the 32-D performance model, the program marched up to 
the beqinning of the transition liner. At each marching atep, 
the proqram performed ai heat-flux balance at the wall and, 
thereby, obtained the wall temperature, 


The input units, required by the Viner-cooling model, are 
S$.T., and since these are the output units of the 3-D performance 
model, data in easily transferred, Profiles for U-velocity, 
temperature, turbulent kinetic eneray, turbulence length seale, 
tuei-mass fraction, total fuel, and CO are required along with 


information about the Clow conditions ino the inner and outer 
annuli. Sinee nearly all the liquid fuel has evaporated at the 
point where the marching process begqing, no dnput fer the droplet 
evaporation ia needed, 


The input data sheets (Miqure 19, Vo sheets) beain with two 


cards devoted to case tdentifteation (additonal tinput informat ton 
ean be found in the input sheet ferma located ino Appendix A), 
Next, 40 cross stream arid cells are speed tied. Tt was unlikely ; 
that this would be suttiedent to obtain a arid-independent solu- 
tion; however, sinee thin wan merely an example case, the maximun 
dimension allowed in the existing deck was setected, Axisyin- 


metric qeometry, selve speeder equations, two- equation  Curbu- 


lence model, reacting species, wall temperature, and enthalpy 
ealculation completed the napecitdeations on Card at. Card vo indi- 
cates that the vardoun protdiles will be printed every 20 mareting 


ntepa. On Card 4 one sees that the macehing process wilt atart at | 
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Figure 19. liner Cooling Model Input Sheet (3 of 3) 
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0.03937 meter and proceed to 0.073 meter, using a step sise of 
0.02 times the grid height. In order to obtain the best estima- 
tion of the radiation flux, it is necessary to analyze the entire 
channel height, which means the inner and outer boundaries are 
walls. XEND and XOUT are, therefore, set to some. large number. 
The variable POWER is used to distort the grid, since more nodes 
are required near the wall whose temperature is being predicted, 
This means, of course, that two runs must be made, one with the 
grid nodes concentrated near the outer wall (POWER <~1.0), and 
one with the nodes concentrated near the inner wall (POWER >1.0). 
If POWER equals 1.0 or -1.0, the resulting grid would be uniform 
in y. 


Two straight walls require only two boundary point specifi- 
cations, handled by Cards 5 and 6. ‘The three radii specified on 
Card 7 are also easily understood. ‘The name list is read next 
with the variables specified explained in the name list input 
sheet. ‘The values on Cards 9 and 10 pertain to the tree bound- 
aries, which are not used on this program; however, entries are 
made for completeness. Cards 11 through 20 specify the initial 
profiles. The 19 points used correspond to the 19 radial nodes 
used by the 3-D performance model and are listed in Table 1. For 
this example the 3-D output was strictly used; however, it is 
usually the practice to combine the 3-D profiles with some addi- 
tional intormation, if available, describing the protiles inside 
the slot Lip. Card 21 lists data concerning the annulus flow 
conditions and dimensions while Card 22 is blank due to the 
absence of liquid fuel. 


The output from the program is shown in Figure 20 (3 
sheets). Sheet 1 begins with the titles, some control indices, 
and the values of omega for the transtormed cross-stream vari- 
able, Information about the annuli is printed next along with 
the initial profiles. Whe printing of the variable arrays starts 
with the value at the inner boundary and continues outward, 


103 


68S 

68S 

aszT 
SOET 
OSseT 
S8tt 
Stet 
09¢T 
STET 
S06 

stot 
OOeT 
SGeT 
OcET 
OLZT 
otet 
SETt 
68S 

68S 


dWa 


0 

0 
LztZ0°0 
87Z0°0 
GTZ0°O 
s6te’o 
99TO°O 
6ZT0a°0 
ozT0°Oo 


€12£00°0 
91600°0 


Tvt0"O 
BS5T0°O 
z8Tto°Oo 
8610°0 
s0zo0°0 
06TG"O 

0 

0 


OON 


L 


0 

0 
€-ACL“9 
€-909°V 
€-ACEE 
€-dty °c? 
€-d89°T 
€-380°T 
¥-a9°9 
€-d98°T 
€-AL8°T 
€-d9T ST 
E-SEL*T 
€-desre 
C-AT9°E€ 
€-gaLT°S 
£-S8°L 

0 

0 


naw 


10) 

0 
zoco°o 
86Z0°0 
v6z0°0 
T6z0°O 
8820°0 
£920°0 
€£70°O 
LTT0°O 
#STO°9 
6£70°0 
s9cG°0 
69Zu°O 
04729°0 
0LZ20°0 
€L£z70°0 

0 

6 


9 


WY Joy aze san{Tea 1070 


S-3729°L 
GS-dZ9°L 
¥-FL8°2 
v-3zZ9°E 
v-9ST°P 
y-Acc"vp 
¥-ALB°R 
¥-ALTS 
v-dzL°s¢ 
P-AbL-Z 
y-SITe 
¥-499°% 
F-dT8 °F 
v-F9L*Y 
¥-AC7P 
F-as9°£ 
v-a98 °C 
$-329°@ 
S-a9¢9°L 


at 


"ZT pue 9T = 
Lt = 

Le°9 L°St 
£z°9 L°St 
6°2T T°92 
€°?¢s 6° SZ 
g°s9 0°€Z 
G°TL c°6T 
B°SL 7° 9T 
£°€8 £°ST 
ect 0° TZ 
£LZ pres 
9eT °° Or 
G°ie 3°12 
€°ed £°02 
L°TL 7° eS 
L°s9 S°SZ 
6°TS c° 82 
c-et T° 92 
£z°9 L°S¥ 
Le°9 L°S? 
au 0 


"TIdWWXA ONITIOOD WANIT WOd SATIAOUd IWILINT 


*{ Piadva 


I jo °Hae ay pue 
*‘£ = ¥Y wzaozZ are sanTepa-n 


o000°0 
¥Sz700°0 
80S500°0 
Z9200°0 
9TOTO°O 
0LZT0°O 
yZSTO0°O 
BLLT0°O 
TENTOG 
98¢2c0°0 
0*S70°0 
¥6L70°0 
8FOe0"O 
ZOEEO"O 
9SSt0°0 
OT8E0°O 
y9070°0 
sit7v0-0 
ZLS90°O 


A 


104 


A oO DOA RrFP WWM tH MN A 
aoe 


ma ON 
aod 


oOo vw 
ed 


rt 


\ 


=) 
eo 


{ 


t 
e 


n Oo 
a oe 


TE te ne ne 


Bern 


"(€ Jo [ 38e4S} 3nd3nO Tapow HutTyTooD) AsxuTT "02 ainbtg 
: 9 $0-3029°2 $0-3029°R 90-399E°2 YO-FIGESE 90O-3L06°E£ 90-30EES? 80-3452 °% 40-3£62°9 49-30EL°9 
95-3L9E Te € 9O-3L16%< «990-3900°F 90-2226 °% 990-3 ebB°S «90-3961°% 96-39T05S 20-2658%% 90-392497% %0-3905°9 
$0-35EE°9 2 80-3290°%R 9O-FESZ°E 90-30ES*E FO-de¥Z°E 99—-3520°E 9890-52999 =0-3590°2 990-3826°1 = 90-59Z0°1 
‘ $0-3029°2 Lo GG-30c9%2 §=§0-3029°L §=$9-9629°2 §=SC—90Z9°L «=SG-3029°2 «=SL-3929°R $0-49029°L $0-3029°% *o 3) 
“OC GGezR0LES 04302279 IOFFELI OL [Oe2es9°E 1003022°S 1OV9e3TOL°F TO49L9T*L 109206294 10037ES%% 
16+3599°8 “ET 20¢5920%e 25e5koe°2 «206907591 DeRELGTL 100S5S9%S [003265°L 100996972 TOSK9CETL TOC 2TSZ 
10+39¢69°9 "2 LOCFTEC Ty ([Oe3Z95°S 1D439LGL°S LOe3EVEE LO*esuo0°Z2 [OesSZ2°] 19¢38EO°L 0043599°R 00+36605% 
09+3022°9 “9 CC430L2%y 09302299 0003042°9S =OOesbze"S = CUCSOL2*S £O30L2°9 00039022°F 00030472°9 °c as 
2 10-90ZE°2 1C-30ZETZ UG-FESa*E WO—9le69°X 10-32G69° I1C-369641I [0-22E6°R 49-3Z215°! ToO-3BEC*T ; 
[5-3269°T Lt iC-S¥s0°t [C-36Ze°T §=rG-31G2°I Id-s609°T YU-3909°1 [O-seIset [C-3929°1 10-3299°T 1o-39g87 3 H 
c-3t29°T LT 10-300S°% G-3915°2 1O-30€S*L 10-305S*% 10-409S°I 16-3€59° FO-3198°1 T0-3ES0°2 O-de22°2 1 
16-302€%2 2 [O-40cE°2 {C-302E*2 TO-302E°2 T0-30ZE"e ie-sOcE°c 10-302E°2 I9-302E°S TWO-30ZE°e 10-3702EC? w29ax0 j 
"9 sc 0 Z0-32€5°% 20-3008°1 £0-3520°2 20-3S7071 =20-3128°1 =20-9099°1T 20-386%51 j 
20-392€°1 *T €G-3E9E'R EG-3400°@ ZO-3Z201°1 20-392°1 cOmsduc*! db-d2as*t 20-3949°1 20-31Z8°1 20-3056" i 
23-3260°2 2 G-3Se8%2 29-3242°2 29-3622°2 20-3546°2 «cO-F2A 22 STH 994DTZ 29-3L04°L = F9-3002°O E5-IGEs*Z | 
t "9 *3 °c “9 *0 *6 °0 *9 9 *0 “0 u3 | 
q 
“9 “9 “0 €0-3692°S EO-3422°9 «=EO-351S°% €O-29LETE EG-35ES°2 EO-FIEO°I £0-3905°T H 
£9-364E°1 SE EC-35907L €0-4999°E 9O-3519°6 30-366c%s SO-5490°R EG-326€°T E0-39TL°T €0-368C°2 ED-1999"¢ 
£5-3C29°2 "EEO stl E ED-F1LGE°S €0-39505% EO-S90e°S cb-dSec°9 ED-99C0°4 €0-9399T°S EO-3829°2 %0-3009°9 
ize] °9 "9 Oo "9 *¢ *G *0 bd 3) *o “9 13n3 
0 L020LS°% IG4sOeS*? 1063069°F TOe3RSEL*2 WO*|BLOL°Z (Oessee°2 1009222°2Z ToO+3gL0°2Z To4+2ZOT°Z 
TleFSEST2 LO*sGZa°E 1O097OsL°o LO9FLSL°Y T96FZSL°2Z Weosssa?i COeseesel LOesZzS°l TOsISes*s TO939eL°T ToeIseerl 
TO03dTL8S 190342252 1Ge4elHee 100I9L1S*e =$009b9s-2 100590S°2 l0+dLe77< Yes 362972 WelE2Z°E (303G6L°E 100300E%% 
$90 59LS°S TO0F0ES*S 10495026°9 10630LS°% 1003025°% TOesbss°s 1003025°S FCes225°S 19¢3025°9 1909945°% *t a Ye} 
20-3225°9 25-AZkSSY 2C-499ETS ZO-99CE*y ZO-3ZTG°E ZO-3EGLTE 2C-4006°E 20-3FDETE 20O-321T*E 20-39726°2 = 
20 299E92 «20-3245692 20-4897 2 26-3092°2 9 ZO-3EUOTS 20-FeeoTi 60-4962 7E CCH 3979°T 20-3219°2 20-3EaE°1 20-309e°T 
2oMSEaTOL 25-3ZEITT EG-s¥ocrG EC-F8S2°S EO-3BSTESL ED—Fe29°9 LO-ZICG*S EO-2h29°S ED-AGTI°Y €0-3059°E EO-32S8°2 
CT-9GTE*2 CO-FE2E°T EG—-3D0S°L FG-3020°T $0-399E°L YOmdcde®s BO—5c4S*C VO-FESTSL $0-3055°2 "0 10-305S°2 Soatta 
‘ “0 €0-30sE°1 “9 “a *s a) °3 9 
i “9 °C *9 “0 *c “6 "O VO-FE9T*9- -1F TST Y ZO-FE7S°T oOaeL T on 3e 
i °g 3eio*e- ba) *0 "C °6 “9 9 
: °G °6 *o *9 *¢ “¢ *O €0-3908°8 effIsIfw t6-3699°2 =Ginva I owin 
; COMMONS L 260355952 CLeaGTI es *9 *0 “6 3a eas) 2 
‘ “c “0 €S-a6192 ZO-39E64Z EO-FLZETT EO$3GGETE PO-Flvd cs EQeST ETE SEF 23NT3 00030994 o1Nl13N 20-24E6°E one 
°O uast °G eGssjed TO-3ZEE*1 =13e *O ejue *o elue 
: TC-se8E°R oatse "O eIS GO-3CLETI «90G Lo onar 1 outa 0999T#1NO} 900DTeCNaT COCOL<x¥1 C e31$! 
; {5-39Z21°E 20037098°S 10*3088°S 3 
12-JETS*Z 2003003°S 10030L9°2 i 
eeoeseece eeeeseeo ss eo eeeerres 
G¥ad NPS viaa Azaunga | 
QO°90D9T OGe50SC°L 1C-40Z2°@ 19-se0772 % 16-3299°9 £0-3€62°9 10-3290°9 10-3999°S 
10-2909°S 1O-3182°S 1G-geCacs [C-39Z21%% lo-seee*s & %-3s252°£ [9-FOGT°E 10-3920°€ 19-3196°7 
1O-9299°2 W-2Atk°2 WC-390S72 (0-3299°2 19-399EECS 16-3460 %% Ube zaeest «15-Isoeet 19-3269°% TO-F6e""7 
L-9902et 2C-F295°s eC-3sCOE°e 20-365E°S 20-30TL°% eU-36EE°R 2G-3976°S «= 6-3909°1 ba) °c v9 300 | 
*c- °9-0302000°E 90*3001°2 YIeJOOETE 2005ET OT | 
utA TIA Te) tadjav ldtied SS3ae : 
ec-2otece "0 2903060°S 200 3060°S we *9 { 
12°Wa @1 ei an m H 
euLssva 2 ©1300" 1 exedol 2 5n2300" { 
- “ge 24 i 
; ISVD Jtawws3 13584 Sad 1985 eb 


ARAMA RATAAIAAL SADA OREASAEAZASASABANSESCASL CAPS CALAIS AALAAAZAMAWARARAMALASARACARANATAR 
BE ARREL MARA BASIRARASALABDADAZAZASALASA SATA MAUL CAAA PASATAM ALIAS SAARLAAMAMALARAMANRATARALLE 


AN now oer nam ney wr 
: Der aou goo icra porary) ood 
oa ere pia eee oor pit 
vache et a sas Mode eae Sted hee ae 
ou wm wen Pwr BNata ves 
Nn Nar eo Oe eo FPO ee 
rn vom TO" oun rar oRN 
ere eee @@e ee? ave ees 
Onn onr nde Dea moto www 
aN NOAM MAEM HAPS UNION ADM AT 
va YO09O 8vG80 BECOS 93000 8030 2BO0n 
4 ore e078 Ltt eee ofee Fite 
rire) es ee wom Peneen el 
nw? SAMO whee unm vee ana 
mae SOAP NARA BAO rugm Qoom 
ow SNOW BEPH FrPQ wwe OMe 
eeue esos esces eee oeee eeee 
ONRD Ane 96 OQ HMM Wi atae ch were 
naw RANA MAG Remy RAMON Kee Lem ‘ 
on Bane Gece cove AOEA MOAA MOOH -—_ 
va eoe¢ee @O000 tt eeroe eee biae foe) 
aes ith betas Ahi HULL BULE Bee 
tn Otero woe GYRE oree were awe 
we AEE ROWA: Roweeyr PEOKR AAP CAVE aed 
a CFOWE GAH OA we AHO CEH BPOM 
eeteant ones eeoee eevee @eaene eeee eevee (a) 
BNA COOK NFO Nene Oe ee er 
N 
wt BAAN MRM Hm NMAN we Cen 
eo E000 Eoas coee ACHA OOS NOOO 
et eee See Fett eee oof Vii uw 
rere wee Us eee Ue me 
SF Svcd 2AVS fOww chad vee feoe Y 
he FHER NOON Cam OREO CUES BHO? 8) 
om ©OFRO RNG PO“ HOT SOS Kine Rect 
eeee aeee e@eet @ete @eeoes oeee eee 
ON RAN DOAN Wen Oe De Me ol (tp) 
_ 
RNY AMMA MED Heme amma woven 
©2090 VOooO C0a0 3000 VOOo ecco ; 
1b @@#@¢ @69% tebe eoeoe trot v 
Pn bbe Md td eu re be} \ 
Chet CONS HOO OKC aco eIadl A 
@IOC COMO AHOO HES eres weary Qu 
VOW OFA Game FoMme wines Wa ty + 
eoenae vee esos sete e¢eaoe eve 
ONNN Ketetet Oto. mem Crane OR mene s 
NAN NOMM MANHW Nee Amem We nw weed \ q 
999 gece 0000 9999 SVE OOH Onoa | j 
itt e@06 +0084 vt oooe #e0e iil ) 7 
BhW Wiha WE Woe WOU Wee HOU wv ‘ 
See CONS VERO Cane ROWE CEE Ce Rm Le} 1 
Ar ER Atm Aon ewes POOr auve ” 
ORR SAME HHKOO FOVE Over vVE WS AA QO 
wees eee ees et C086 oovee eee oeoe x 
ORAR Rene Detet NAM Came cet SARA : 
AWA NEM mean eee nem ee eae oy i) 
99 coco ooce Coco ecee aaac cane Cc 1 
‘ ooo Soe iit ooee soto Fait 4 
WU Ue i le MUU wee LRUe ae ‘A 
OO Cetin CHO Orad VA me CUA Pam ed 
Srl SSO NON ogee VOCS MUR mame ‘ 
OBR ORVN Ben FON  Suee vas Rane 8 : 
weno see shee sere ore e e8 be Stee 
OMe Cet Cmte LAP Cmte te ma si) J 
q 
i 
ANA NMAM MOEA Know Ree Orn wee 1 
O08 800 S600 cose BsnEe KAM AOOlO My 
144 @e¢0 @606 ttl ooo eee bliit y : 
Woe WwEbw htt $3 to VO Om a) ee ae x 
nee OW TROO Gams NAN ww APIS oy r 
NAF BAOHVOH NAN PF Mo HREH SPPFF eng s "4 4 
ROW Bet Gan FAAS WA AME. Aa! 
eee ante even see8 eeee sere wees ro] 
ONAN Rete Oat omer Daieeet adewetee Pomme 
NAN RNR MEH APES MAMM Swen wwe 
Vea 82900 OO090 vouD OVO VOWS AAVE . ‘ 
bli ee weet FTtit eree oF oo THIS Qo 
Www Wl whe lel eer ren a emo eh] 
NWO COMA FANS EMOR MOwe Pe miter ™N 
we PORN NINN HOnM FOP FOO Oey 
OOD FEN DOME Care WOW WMO. OJAw 
@eee 6888 06a nee eorrne e800 e844 y 
OMA Avett OAs Nah Roem tees Fy wy 
3 
NAN NAMM Men An ST RNMM AWM Awe y on 
eo 8EV0 CYV00 ANOS BVOO YVoow Sooo 
4 @eooe Oo ttl eoooe 600 Feil -_ 
bbe S— 4d hp- dill b~eh, d-glllibekerterderdl Dhak hdl tende-dhadiiind Load fa 
seo OM VNON ONT Omer AMEN OAM 
$$ Fro NOFOD RHOD AOE POSP NLM 
© ONE O00 FLOM Ova ww Neo 
@eese eee eee eee @etee ees eee 
ORNN Kew ONAN ROSH One tet PL ifiw 
MANA NNMM Beno anys ANIA MAW AA 
8o S900 ca0cOo soce Yeo S800 YvuovU 
Cb ¢60% e068 131 eee e066 Fils 
Ce tL eT lc Wiis We 
nan OWLS Sor man HAN CANN FLO" 
Ren OeoN AAHO CF Ay Hate NOSE RO e 
Pen SAMM nor CMe Vevey FNOHH Ooo 
eeee @enee eeope a@eaes e@eesn easee eee 
AMNN AO QR Mee Siem Ate SRE 
- ~ 
-_ -— 
a - t rd ’ 
- = a 2 Q 2 } 
zr wi a w a £ 
rN ~ « a id a a 


106 


oad mon niga) ray tin lenand ons wry mtn Nm Doon avt Dal 
oO ano ooo 096 eaa ea0 ooo eco oao0 ooo ooo 6e0 oO 
‘ oot e+ poo rts pie ooe toe toa +++ tee bad ‘ 
wa dd ad ‘as Las ee Pe aed ww wid Wh Sa end Ae Wet hs ad ed hd tase wa 
w one Pear aor aan 900 ere em met ant Naw tn u 
+ oma BHO Ram Te cro Par ove moo Rao Dow om oe 
1 o end am or mp4 une re TAS WY NO aom Yow mew ea 
t oe eos ene eee ese eee eee ene eee ese eee ave ee 
on Ne um Non onm Netet NN mano ANN Dott atom oto on 
at 
1 Dal wren mn 
1 Q Qoco ° 
: * Piad ° 
bed wl wu 
n CAay Hoa 
NN S Ah wou 
-” Pmimed 
ee aeee ee 
om Wr oe ON =» 
9 
mo awn moO uw 
a OO eneaog « oO rt 
Oo ts #ttt @ t@ 4 . 
1 Wee S ntty-4- ee 4, ~ 
w Ov TNO w OF bes 
st A Re CEme N FRORE isa] 
at - me Cf6N mame | 
ia m ee etea mM 8 
ae oRN were 8 ed 
a a = ie] 
oy is MONA arma MEPS OMY mei Om er ONAN RMN ARERR mem ~w 
mT eo 2066 S000 SONG COCO ooo GCooo 9F29° SHOR 98000 F000 ra wn ea foal 
| w Fatt pose tere vile tht @¢s¢4 Ter «ttt #000 eoee Fei wt 
=e erraeel emmy) ave Cem CeCVeL OMmETCUOLVO] Cam CIUQT Cle mma] vel Cav egmmrelvey) Clvemmmrrl env eive mmmvel very vemmmVaveT OY OMmmYeD cre] Tan rey 
ve ES shee ARID FEE FIP SGCNM CEM AMO NN SARI EO AINAN ON YO 8 ee ww 
B= BOOO OEON FORO MNO Dee NOVO VOM O8em RVR CAN GOED mOMM &@ mm GOL 
¢ VeON OOO CONCH MEAL MNRM CRAY Mem MONO OWMHKD OOSSD ring ra Vv 
ee eene aoa esos eere eaene e@eve ever aecoe eeeoe eeoe eee es 
NO MONS NOE MME NESE Ne ON BOOM MINER Ree Need MMHG No 2 
zs PANN cident PPO RB IMRD tee Deis RAP RAI RI AA Rem oa 2p) 
oO ©2066 ScEeeeo 9808008 oCA0 8000 2060 Ge6eo eeoo Geoo oogo a ~ 
wl rt eee Freee tees eet e+oeoe had tt b ooo Sees Fine wl 
eta Wb ls Ul el sok 
OIF Fees REHO MORE BOOM CRON Ce@e OQietwm NEOM AVOT NEOCO ANNAN COMO OI oe pa) 
uO C9060 WMOO COMMON aE FENNR NOt BOP CMR TvOm NeOt DOGe FY MNNO Ww OCC 
tes MPFO eWeEM GWMM FRM MNO Ot'dn) Med VEU SLING wot Twos ao 3 
eceve eaest eevee eave eane aere eeee @eae ces eee eee e@enee awoee Q 
mF PO Mme ENR I Mees Cute ROD DOM BRUNI eteiet mint Piet mh a) 
Lied | aed 
we wre io} 
aim CONN aetrm REEL OOMM amet Ome Wet CANN AMM ANON Meme alm 
° 9995 9999 ACG FOG0 SOS 9OGEO SeGe S009 GOGO Joes PgHC c) io) 
eat Tlet #40 45 1 thy tedbs wove thet PEDVE €440 eee ot 
ecu ee oe) en eee eo ae ow 
ImO © ¢ ee OC ORR err CR HR CHAT YEO OE BOOS HUNK ANA CK wO see ed 
ube OCC SLCH QRKOW YORE AHR Nemes Avr COCK FESR CORE HET NAKRY wet COCs ri) 
cam Cah MEMO Ref HOOD MAE BACH CHEN AEA BARE GOOD Nerd" at 
mo oa. ese e e8ne bees reee F846 hte £888 FPO FCAe8 sees Fee oe Le] 
see CLAM PRU TUR ee Keedet Mehmet BOE ONAN Ome mito oP FU: s*e e) 
« o 
me oI 
VORA mete ALIN G RAO mtmtetedt miotet PPP RCIA QUE nh ee ue = 
e O20 ANCHO GENS GOGE SOG ANCH COOHE ANGE MOSES een GooR ko 
eal Pert eeods Feed et dae Fett aoesn Fret et st Coat +640 FLEt bet on 
weuw We a a a lk Ul 
COM 2 eee OHAO KAMO Lael HOOD COSHH HLeY PROF CHEM Cher RAIN FORO Com sees = 
AG COCH Meetin -AA FRRIN MOAR NOLO MVOR FRE Be Se CHM. CcCHeEer COA af oo¢r ved 
nic CePA C8EF OLCTE Comm MEM KOOP OYNN MEMO OWE COS Ohm mi 
oe eee vees cco e aeee sete Oho Bees HF ee Bere 2488 BO08 oe Com 
enor SEEM HOUR MD et wiRintet Alte mete BPR IA AIR! Downe mteteae Mem 8 eee 
we ro 
we ur 
vem VR NA etm FEE OEM meee mBoem AEP RNIN fUmmr INA Yi Wm 'o) 
wie onee HOGG SOOO GOGO GEG OOCOF O00 AOOCO 92OH06 HOSO Gone wo 
m~oe orb ooo Fret ober tre eooo Feet Phd &eooe e444 FREE mae 
au MW. Win Wee Ul MWMUW Ut Wh ek aw bet 
© ™% FS eee ARNE BNE YRMN OWOM Ow SOOT NOSN YUBO BMMAW AAI NWO § O sene 
29 9000 BON ALAS FMM fem Nee NOD MMOD BWOW MRO Gewow MEN Zane COGS (0) 
wow SNA BHO OPES OOF MATT AGEN ARLSO CMAQ LMR 1000 NNOR wom al 
wioe eeee eee eanoe eeene aeee @eue #eoe eeoee aeoee seen eves miee 
anon Om AWA NNN oP Pe Art ttt MON PR NINA Pitts ntetoton mmON wNON | 
Omi Oo 6 Ont 8 O wl 
om ' 4 om ' ‘ 
QOmea wou YOUN Ae OM POE mini tei PPT ONIN AMM We ee OMB uw 
Qem Pw 9NOG% 2600 ONG SEOVO G060 99000 Oooo soo ©0000 G900 19600 See F 
“me Crthl eee Lead Ford Pett wees tte ued Pose eee FEL dt meme nt me 
mt oY Aka ll ll ll ry bi 
~~ eee mMNAMm FOC GAO FORN CRPeN fener DANS BOG NRAN ANRN OU YS ™ @eve 
DeD mOFO TN ASEM HRN OPO ARO Cmimm OPO NTN ROMY OMoH ama es moc Oo 
Omw t OK OLD FEO wOhe SEMM COMM Nea MmEOQ aaVWe GOOG Cdden ves | N 
ke @eee oF oe 888 hee one CHF CHF ete Cosh eee FOO e me 
a . 1 AFM NIN OMe NAO Niet SRR EMER NINN Det cites mermOUN o eo 
oof # ~ a cv) 
oou7 Ss eoaosn7 
o*+7 = JAM Aaty i OMIA POSS ete Ome ALP ON NNN nmin Aer GO & ~ + M4 
oO wa Oo So00G COO =OGQ 29690 90013 O00 YOCOD SY00 BY0O VOOU S000 Ou © :: 3 
tet Pet oOo FEE H FEED FHL Shoe FEED FL tt oot +oee PEt aowrie ws 
oem 9 ui WU oO 
ana a rorn MO NAMM ANFRM ONE EC NOMS SIM CHE Dede NNN LOND 7 © a wt wh 
z ° OTH Oe WN OOO NOs Meher AYN ATCO SED 2HOH BAMOD F£ e ae 
wer Nm were Qe NAN BOF MOME TOR BROW mar A NOMS OOOO MON Wet Nm fy 
=o oo eee ©8 se 6808 e©8ee seve sees Bees sete snee Fete Bete mO o °o 
hos WME ARAN ARM PP NIKI ION MO Rt OM ete crt Sw) to 
oO fA ow & o 4 wow 
Oaw oO + om ao + 
Oud Pow manuel Metet OMAR OME ninieted ota IPT BRAN je mA AAP OWS Ft 
OL mm oO fellate] QO 900 2999 soo voo ooo 9000 GBS9O |8GOU SOS Sofa —~ wb 
mu 88 pied toe titbt ebat tad ooo tet orhed oer oe (tlt awe 2 8 
sae & Mas as Deed wall <dvads hash hpntredy <i wed= deeded iene] Aaya ks Ale Wie el 
m omaw See SOFos mM Sora ron ator mane moO NNN NOTE mM 
« rare hast MAO OF ON Nae uvo Fae NOHO Sts uDwOM LarN -« fw ‘ 
Loa] L} LJ wane NTO OF On FON Mat om 2020 oom SAF OOWHD —OMM wi O 4A . ! . 
iyo o eeee 2000 885 Bptee eeee Cee 48988 wee eee seee seee ' ace sj 
- NAN OF NN ete atret fuNee ON Oe Runt tet Dedede FOR ml ew | 
> 3 2o 3. 2 j 
aa o aa 
To ed : 
Oe i 
‘ yt om wn ra - a ~~ et -- . 
. w ool a 
a » 6 bad = w a - we os 6 
a] = om . i » =z -_ & o J om st Fm 4 
= o = o - w w zr iw) < Tc wer wow 
ua & a > is pw i=} a a - a a - w ‘ 


107 


Since y is always referenced from the inner houndary, y(1) 
is always zero, so the radius of the inner boundary is printed in 
its place, It car also he seen that the y values are more closely 
spaced near the inner boundary, since this was the wall tempera- 
ture that was calculated during this run. Sheets 2 and 3 of 
Figure 20 show a printout after 80 marching steps as indicated hy 
ISTEP. The program has traversed to an X-location (XU) of 0.0706 
meter, which is very near the end of the combustor. At this posi- 
tion the inner-wall temperature is 958°K. The other profiles are 
also printed with KE being the turbulent kinetic energy, LEN the 
turbulence length scale, PHI the total fuel, RAD the radiation- 
composite flux, and AMU(I) the effective viscosity. It should be 
noted that the boundary values of the radiation-composite flux 
are never used and are therefore left at their initial values. 


TRANSITION MIXING MODEL : ; 


The 3-D performance model was used to peedict the flow field 
up to the end of the combustor liner (Figure 15), but it is not 
capable of handling the geometric configuration of the transition 
liner, For that, the transition-mixing model (TMM) is used, 
Using initial profiles as predicted hy the 3-D program, the TMM 
marches through the transition liner, thereby predicting the tem- 
perature distribution at what would be the turbine~stator inlet, 


An enlarged drawing of the transition liner is shown in 
Figure 21. The inlet and exit planes are marked alona with 
several intermediate ones. The location of these is a matter of 
user choice, but should be enough to simulate the actual liner- 
wall geometry. The radius of curvature is constant, In this 
case it is 0.5 inch for the inner wall and 1.9 inches for the 
outer wall. The input sheets are shown in Figure 22 (3 sheets). 
(Additional input information can he found in the input sheet 
forms located in Appendix A.) Card 1 is allocated to case iden- 
tification, Card 2 shows that 40 cross-stream intervals are se- 
lected along with axisymmetric geometry, K-¢ viscosity model, and 


108 


6.0 


ee eee =, 


Transition Mixing Example Geometry. 


3NV id 1FINI 


INV 1d Lixd 


(E) 


2.0~«~«T*N 
ee ee eee ee Le eI 
F 3. 
(X~ IN) 
Figure 21. 
109 


CASE TITLE CARD (12A6) 


PTRANSIITion Mixing |Mobded | Exampde a ae 


CASE TITLE CARD (12A6 


CIT Ee eC eras 


(8 (12, 8X)) 
KRAD MASSTR ISWRL MODEL INERT IFLUX ITEMP 


me ler ee ee ee Ter ies er 
NSTAT NPROF NPLOT ITEST LASTEP 5(%5, 5X 
#1999 Joooas]ia1971 [3109 LTT 


(8E10.4) 
2ur ZUE ZULAST FRA ZEND ZOUT PRESS POWER 


Lo7see Jovsec}23e J.o2 [1%  jio.  frog2S [io | 


NBP ICURV NRCVI NRCVE (8(12, 8X)) 


eee (ae (Se a ee ee 


(8E10.4) 


tr 


w 


> w 


wa 
ec ata te te ee dae 


(8E10.4) 


aI) RCVI 21 RCVI 2i RCVI aI RCVI,... : 
5 Saar ae RE oP RET ARIE RET : 


Figure 22, Transition Mixing Model Input Sheet (1 of 3) 
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(GF10.4) 


{ aE) RCVE, 2E4 RCVE, 2E, RCVE., ZE4 RCVE,... 


(8E10.4) 


RA RB RD 
a .2s4 fare | 


NAMELIST 


9 |ginpur ]icreendfe | TT Ts | 4 


UA UD (8E10.4) 


lo [ee Cen ee Pa rane 


SS Ean SRR Se eee Ren, eae Pace De 


NIN 2 
oa ks eee Ce es Ree Cee 


¥ VALUES (8E10.4) 


\ 
| 
| 
F2A F2D =TA to Twall (8E10.4) | 
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l 
| 
! 
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U VALUES BE10.4) 
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TEMPERATURE VALUES BE10.4 . 
16 | | | | 


Figure 22, Transition Mixing Model Input Sheet (2 of 3) 
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Figure 22, Transition Mixing Model Input Sheet (3 ef 4) | 
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nonisothermal flow. Card 3 will allow a printout of the variable 
profiles at every YO marching steps. The roitial and tinal 
values of 2 and the marching direction are placed on Card 4 along 
with the control on step size. Since the boundaries are always 
walls, ZOUT and ZEND are set to a large number while POWER equal 
to 1.0 torcen the initial grid to be untform in y. POWER »1.0 
will distribute more nodes near the inner wall while 0 < POWER 
<1.0 will distort the qrid toward the outer wall. Next, the num- 
ber of boundary points and number of raditus-of-curvature pointa 
are specified separated by a Flaq to indicate that radius-of- 
curvature effects are to be included. Card set 6 and 7 reads the 
Actual boundary and curvature values which are listed in Table 2, 
The values listed are in inches and had to be converted to meters 
since the TMM requires S,T. units. The initial and final values 
of radius-of-curvature were just chosen to be a large (1000m) 
number, Cards 8 and 9 are quite obvious with the name-list vari- 
ables listed with the input forms. Cards 10 and 11 deal with 
“free" boundaries, which are not used for this program, so the 
values Listed are only for completeness, Card 12 indicates that 
19 points on the (nput initial profile were used and correspond 
to the 19 radial nodes employed by the 3-D program, The various 
profiles are read on Cards 13 to 21 and are listed in Table 3, 
These protiles are merely the exit plane profites for #0 = 6 
degrees (in line with the primary orifices) as obtained from the 
J-N output, 


The output of the TMMs in tllustrated in Fiqure 27) (4 
sheets). It begins with a list of control indicer and tmportant 
quantities followed by w, the transformed crosn-atream variable, 
The specified boundary valuen of X and R for the inner and outer 
wall along with the value of 2, the dintance along the wall, are 
linted next, Following these are the values of radius-of- 
curvature and the initial profiles, Printing of the variable 
arrays starts with the value at the inner boundary and continues 


eutward, Since y ist alwayn referenced from the inner boundary, 


cand as cen eer amine ai — sae 


ae ee 


ae ee 


TABLE 2. 
Point XI 
1 0.0 
2 2.9 
3 3.01 
4 3.10 
5 3.26 
6 3.37 
7 3.40 
8 3.37 
9 3.26 
Lo 3.10 
ll 2.90 
12 2.60 
Point aI 
1 0.0 
2 2.9 
3 3.01 
4 4.27 
5 4.47 
6 4.77 


TRANSITION MIXING MODEL GEOMETRY INPUT, 


RCVI 


39370.0 
39370.0 
0.5 
0.5 
39370.0 
39370.0 
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XE RE 

0.0 11.8 
2.9 11.8 
3.37 11.73 
3.80 11.56 
4.45 10.98 
4.76 10.20 
4.75 9.50 
4.45 8.79 
4.02 8.38 
3.51 8.10 
2.90 8.00 
2.60 8.00 
ZE RCVE 
0.0 39370.0 
2.9 39370.0 
3.37 1.9 
8.19 1.9 
8.81 39370.0 
9.11 39370,.0 


TABLE 3. INITIAL PROFILES (SI UNIT'S) FOR TRANSITION MIXING 


| 
| 
E 


EXAMPLE 
Node r u 1m Temp. 
‘19 0.04572 32.6 15.6 5.86E-4 879 
18 0.04318 32.6 15.6 5.86E=4 879 
17 0.04064 22.8 18.0 6 .03E~4 1310 
16 0.03820 26.3 27.0 6.16E=4 1470 
15 0.03556 29.2 30.3 6.50E~4 1530 | 
: 14 0.03302 30.7 30.6 6.64E=4 1560 | 
13 0.03048 31.7 29.6 6.65E~4 1570 
12 0.02794 33.0 28.7 6.59E-4 1560 
11 0.02540 $4.6 29.9 6.42E=4 1530 | 
| 
10 0.02286 37.0 61.8 4.65E~4 1470 : 
9 0.02032 34.5 42.9 G.24E=4 1550 
8 0.01778 30.6 27.0 6.58E-4 1660 
7 0.01524 28.7 26.4 6.57E=4 1700 
6 0.01270 27.4 27.8 6.59E~4 1/20 ! 
5 0.01016 25.9 28.3 6.49E-4 1690 ; 
4 0.00762 23.6 25.3 6.21E=4 1570 
3 0.00508 21.6 20.4 6.10E=4 1370 
2 0.00254 33.5 16.4 5. 88E~4 879 | 
| i 0.00000 33.5 16.4 5.885 -4 879 = 
: K=7 
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Y(1) is always zero; therefore, the radius of the inner boundary 
is printed instead. Note that the pressure variation across the 
grid (PRESS) is essentially zero due te the large values of 
radius-of-curvature (RCVI, RCVE) specified at the inlet plane. 
Sheets 2 and 3 of Figure 23 show the printout after 160 marching 
steps, Along the outer boundary, the program has marched to 
0.1639 meter (ZUE), but the inner boundary is only to 0,09497 
meter (2UY). At this point, there exists a considerable radial 
pressure gradient with the inner wall at 744 N/m? lower pressure 
than the outer. In addition, the velocity profile is distorted 
so that the maximum is very near the inner wall, Sheets 3 and 4 
of Figure 23 show the output at the exit plane. The radial- 
pressure gradient disappears as the radius-of-curvature effects 
are no longer present. A plot of velocity and temperature is 
alao made. The top line of the plot corresponds to the inner 
boundary while the bottom line Is the outer. However, it must be 
remembered that the inner boundary here is located at what would 


be the stator-blade tip and the outer at what would be the blade 
root. 


EMISSION MODEL 


The calculation domain that would be used by the emission 
model to predict the emission output of the example combustor is 
illustrated in Figure 24. Starting with the initial profiles 
near the dome, the program marches using the recirculation zone 
as the inner boundary. The mass and specie concentrations in the 
recirculation zone can be estimated or calculated by several 
methods, and this provides an entrainment rate and boundary con- 
ditions, Once the centerline of the combustor is reached, the 
inner boundary switches to an axis of symmetry. Cooling slots 
and radial-injection points can be handled without stopping the 
marching process. To input the information necessary to perform 
the above calculation would be quite lengthy and illustrates 
nothing about the emission program. Usually, the boundary along 
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the recirculation zone must be handled in a manner that requires 
internal modifications to the program similar to the Cards MA.528 
through MA.743 that are unique to each case analyzed. Therefore, 
for the example case (Figure 25) the program is started at the 
downstream edge of the cooling-slot lip and marched to the exit 
of the liner. 


As illustrated in Figure 25, card set 1 is available for 
case identification. Card 2 specifies 30 cross-stream intervals, 
axisymmetric geometry, species, two-equation turbulence model, 
reaction and nonisothermal conditions, while Card 3 produces a 
profile printout at every 40 marching steps. The initial and 
final x locations are read on Card 4 and have units of meters, 
since all the units of the emission model are S.I. The step size 
was chosen to be small as is necessary with the l6-step kinetic 
scheme. Since the boundaries are always walls, XOUT AND XEND are 
set to a large number. A value of POWER equal to 1.0 provides an 
initial grid that is uniform in y. POWER »>1.0 places more nodes 
near the inner boundary while POWER <=-1.0 places more nodes near 
the outer boundary. Since both walls are straight, only two 
boundary points are needed, Cards 5 and 6. Cards 7 and 8 define 
the initial grid, inner and outer radii, and the name-list quan- 
ities. Cards 9 and 10 deal with quantities at free boundaries, 
which are not used, so the data is included only for complete- 
ness. The initial profiles are the same as those used in the 
wall-cooling model example. Nineteen points are used and corre- 
spond to the 19 radial nodes used in the 3-D performance model, 
Three additional profiles, listed in Table 4, are needed for the 
16-step-reaction scheme. Since for this case there are no cool- 
ing slots or radial injections, the remaining Cards, 24 through 
51, are omitted. 
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CASE TITLE CARD (12A6) 
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CASE TITLE CARD (12A6) 
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3 [29999 Jooovof99999] drag tT TT 


XU XULAST FRA XEND xXouTr PRESS POWER 8E10,4) 
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1) eee eee REA PAGING [Pa (NR a EE 
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Figure 25. missions Model Input Sheet (1 of 4) 
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NAMELIST 


8 [sinPuT | KREAD=|!, PREXIb/ = 2-PEH/C, ARGow/=S#oco, | s _ 


CRIz3.0, CK a10, ee Ene EFus.§, Eo2el0 


VTD 8E10.4) 
9 SAE ENE (Tan FS CT A DE RENAE 


Twall (8E10.4) 


Lo (A! A (7: SC; C7 Ml 


NIN NUL NUE NVI NVE (5 (12, 8X) ) : 
Ts a CC (ee CC ee ee ee 


Y VALUES (8E10.4) 
12 

U VALUES (8E10.4) 
13 

Va VALUES (8BE10.4) 
14 

TEMPERATURE VALUES (8E10.4) 
15 

MFU VALUES (BE10.4) 
16 

MCO. VALUES (@E10.4) 


2 


1, Ee en tee iaie Ga ee Caen Cee 
Figure 25. Emissions Model Input Sheet (2 of 4) 
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MCO VALUES (8E10.4) 
18 
MOX VALUES (BE10.4) 
19 
MH5O VALUES (GE10.4) 
1 Lemar Mere (ae ee (ER eae] [Eeaesieaa [Eee Rae 
MH. VALUES (BE10.4) ; 
21 ee ee a | 
\ 
TURBULENT KINETIC ENERGY (8E1LO.4) 
22 ! 
TURBULENCE LENGTH SCALE (8E10.4) 


hi Ae a ee aS a ae ae : 


SKIP FOLLOWING CARD SET IF NUL = 0 q_eeemeeeee 
X = LOC. OF INTERNAL COOLING SLOTS (BE10.4) 
fee SR Se! SAA ET STE NEES | 


LIP LENGTH OF INTERNAL COOLING SLOTS (8E10, 4) 

25 
U - VELOCITY OF INTERNAL COOLING SLOTS (QE10.4) 

26 
VT - VELOCITY OF INTERNAL COOLING SLOTS (BE10.4) 


a) anne LE (a) (a aR (reine AE) at 


Figure 25, Emissions Model Input Sheet (3 of 4) 


125 1 4 


ot ia os 


\ TEMPERATURE OF INTERNAL COOLING SLOTS (@E10.4) 
q 
i 
FLOW RATE OF INTERNAL COOLING SLOTS (8E10.4) 4 
29 
Y 4 
i 
SLOT HEIGHT OF INTERNAL COOLING SLOTS (BE10.4) | 
30 : 
— 
| A 
bq 
SLOT TO METERING AREA RATIO FOR INT SLOTS (8E10.4) r : 
31 ; 
SKIP FOLLOWING CARD SET IF NUE 2 () cnemmmemcnmenmnnme ; 
X-LOC OF EXTERNAL COOLING SLOTS (8B10.4) | 
32 4 
y 
1 rf 
ie 
LIP LENGTH OF EXTERNAL COOLING SLOTS (8E10.4) oo 
33 
U - VELOCITY OF EXTERNAL COOLING SLOTS (BE1LO.4) 
34 
V,, - VELOCITY OF EXTERNAL COOLING SLOTS (8E10.4) 
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Figure 25. Emissions Model Input Sheet (4 of 4) 
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TABLE 4. ADDITIONAL PROFILES FOR EMISSION MODEL. 


Node 


19° 
18 
17 
16 
15 
14 
13 
12 
11 
10 


9 
8 
7 
6 
5 
4 
3 
2 
1 


M 
co, 


—oentemene = 


0 

0 
0.0318 
0.0368 
0.0311 
0.0484 
0.0534 
0.0494 
0.0274 
0.0199 
0.0527 
0.0657 
0.0596 
0.0537 
0.0486 
0.0438 
0.0400 

0 

0 
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0.0243 
0.0272 
0.0291 
0.0304 
0.0309 
0.0282 
0.0169 
0.0123 
0.0262 
0.0339 
0.0338 
0.0333 
0.0325 
0.0314 
0.0293 

0 

0 


ww 


eer 


The output of the emission model is shown in Figures 26 and 
27. The initial printout consists of some control parameters, 
the value omega (the transformed cross~stream variable), and the 
initial profiles, The printing of the variable arrays starts 
with the value at the inner boundary and continues outward. 
Since y is always referenced from the inner boundary, y(1) is 
always zero; therefore, the radius of the inner boundary is 
printed instead, Figure 27 shows the profiles at the exit of the 
combustor, XU »« 0.073 meter. The profiles of NO show small 
values (much less, however, than one would expect as the prin- 
cipal area of NO formation) occur upatream of where the emission 
model is started. 


FUEL - INSERTION MODEL 


The fuel=-insertion model could have been used to predict the 
droplet trajectories of the nozzle in the 3-D performance model 
example case, Figure 28; however, since that program contains its 
own spray model, a asaimple illustrative example was selected 
instead. A two-dimensional grid 5.0 X 1.5 inches, with the spray | 
originating in the lower left-hand corner and processing a non- 
uniform flow field, was analyzed. 


TO te en cca = = Lind 


Card 1 (Figure 28, Sheet 1) allows for case identification. Pe 
Card 2 specifies the atomizer type, in this case an airblast noz- 
zle with a 10 flow no. and 90-degree cone angle. The other fields 
are left blank as they pertain to dual orifice and/or simplex 
nozzles, Cards 3 and 4 are omitted as they are used to input the 
SP versus fuel-flow curve for the secondaries of a dual-orifice 
nozzle. Card 5 specifies some dimensions of the airblast nozzle, 
the filming diameter, and exit-flow area along with the airflow 
rate and temperature, The fuel type, flow rate, and temperature 
are given on Card 6, plus the flag to specify a nonuniform 2-D 
field, Other quantities on Card 6 are required only for other 
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l TITLE 80 


Ful INSERTION MODEL EXAMPLE 


AIR 
PRIM SEC. ASSIST 
® ** FLOW FLOW CONE SHROUD PRIM SEC 
ATOM AIR NO. NO. ANGLE EFF ORIFICEORIFICE 
TYPE ASSIST(JP4) (JP4) DEG AREA DIA, DIA, 
PeH/VPSI 1n* ss IN IN 
2 


*ATOMIZER TYPE: 00001 = SIMPLEX 
00002 = DUAL ORF 
00003 = AIR BLAST 
**AIR ASSIST: 00001 = NO ASSIST 
00002 = WITH ASSIST 


"FOR ATOM TYPE = 00002 (DUAL ORIF) ONLY (LEAVE OUT FOR OTHERS) 
INPUT SECONDARY FLOW SCHEDULE 


Wg ™ SECONDARY FUEL FLOW, LB/HR 
AP," APsec, onrr * “Prrow piviper vans’ P/P 
CRACK 
POINT 
FLOW 
W W W W w 
$2 
eee. See | Gaal aca \ Geko 7 
| Se ee Se a ee Tees 
CRACK 
hr aa ey OP, AP 
to ee a ay at 


7 anneal Sikes) es ie ae Re 


*FOR ATOM TYPE # 00003 (AIRBLAST) ONLY (LEAVE OUT FOR OTHERS) 


FILMING NOZZLE FLOW AIR 
DIA AIRFLOW ARER TEMP 
IN 


f 
; 


ew 
w 
tr 
Ls) 
we 
us 
wt 


Figure 28. Fuel Insertion Mode) Input (1 of 3) ’ 
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5 ey 
eels Shot re ne Pn te 
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ween AIR AIR 
*** AIR FUEL FUEL FUEL ASSIST ASSIST 
FUEL FLOW TEMP FLOW AP AP TEMP 
TYPEOPTION °R  LB/HR PSI PSI oR 
1 


TYPE | 00004=IP4 FLOW 
OPTION 


***PUEL | Q0002mIP5 ****AIR 00001"UNIFORM GAS STREAM 
00002=2=D FIELD OPTION 


9 
angi Yona, Sanas)  aky; MA ee 
1 ll PPS 21PSTA 31 IN. 4 °R 51 


7 2 2 fee? [SoS 
t 


t t t t 
¥ 


Xyoz ¢ Ynoz ¢ Peas “max ‘max 
IN. IN.  PSIA IN. IN. 


CARDS 8 THROUGH 13 SKIPPED IF AIR FLOW OPTION = 00001 


IN = NO, OF X-DIR POINTS IN 2-D FIELD 


2(12, 8X JN = MO. OF Y= DIR POINTS 


IN JN 
ge Cea (a CP (ARR! CEE (RP RG 


X VALUES (8E10.4) X = LOCATIONS OF GRID POINTS (FT) 


¢ [O, f. 2083 |. 4/6 


Y VALUES (8F10.4) Y + LOCATIONS OF GRID POINTS (¥T) 


ro [O. poéas[i25 fo | en eae 


READ ALONG + X LINES STARTING WITH 
U VALUES (8E10.4) SMALLEST Y VALUE (FT/SEC) 


11 (30. [7s tien fT panerren! Iie 
3S. Ro, lo, 


ca a ee ee a a a 


Figure 28, Fuel Insertion Model Input (2 of 3) 


tie Z-D FIELD OPTION 


READ ALONG + X LINES STARTING WITH 
V VALUES (8E10.4) SMALLEST Y VALUE (FT/SEC) 


2 Sven [= ec) (An (SARS CReme! RREO 


Oo. Oo. (or 


READ ALONG + X LINES STARTING WITH 
T VALUES (8E£10.4) SMALLEST Y VALUE (DEG. R) 


13 loeo [loge TT 


‘ loeo loso O60 


oe |loéo Jtoco | TC 


Figure 28. Fuel Insertion Model Input (3 of 3) 
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nozzle types. The values of Card 7 can have two meanings, 
depending on the type of flow field. In this case, they give the 
X¥ and y location of the origin of the spray, the pressure, and 
limits of the two-dimensional grid. Cards 8 through 13 are 
required in this example, and it can be seen that only a 3 X 3 
grid was used which requires 9 values of U, V, and temperature. 


The program output is illustrated in Figures 29 and 30. The 
first items printed are some of the input quantities, along with 
the specific gravity and viscosity of the specified fuel. The 
fuel AP and velocity are given next, along with the calculated 
value of SMD. For each of the five droplets, the locations are 
then given at the point that selected fractions of the fuel that 
had evaporated, In this case, droplets 1 and 2 evaporated within 
the specified boundaries. Fiyure 30 shows the trajectory output 
for each of the five droplets, The output consists of a pair of 
lines, the first giving the X and ¥ location of the droplet and 
the second giving the diameter, temperature, velocity, and frac- 
tion evaporated. 
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APPENDIX A 
INPUT SHEETS 


ANNULUS LOSS MODEL INPUT SHEET 


1 TITLE ~ MUST HAVE FOR EACH CASE 


ANNULUS INLET FLOW CARD 
NEL NELI 
1 345678 12 M2 9, PTL 4, TTL 4, BETA, DP/P 6) 2 


DOME INLET FLOW CARD (REQUIRED FOR INTERNAL LINER FLOW ONLY) 


NELI 
12678 11" a3 PT 3, TT 4 BETA 5) 


3 fod SS —— 


CONSTANTS CARD 
1 i FRIC, , BLKF 3, TANS 4, CDB 5) 5K 6189 71 TREAX 


1 Title ~ Run ident appears on printed output and plot 


=# Ll some (or all) holes fixed, inlet flow fixed, 
iterate to get preasure drop. 


= 2 gome (or all) holes fixed, pressure drop fixed 
iterate to gat inlet flow (input Wl is firat guess). 


j2 O no plot 
EPIC in 1 plot drawn 


= © output printed after converged solution 
IDBUG : 1 output printed after each iteration] Use only to de-~ 
= 2 output printed after each element | bugg 
non-convergaence 


2 3 NEL = total number of element stations (card 2. only) 
NELI «= element ID no. (NELEM) at annulus or dome inlet 2. and 


ITER 


3 
Wy and W = air flow at inlet stationa, lb/sac 
a and PT = inlet total pressure, PSIA 


and TT = inlet total temperature, R 
neta, Beta = swirl angle at inlet _ 4 


DP/P = total pressure drop, PSIA/PSIA (first guess if ITER = 1) H 
card 2. only os 


4 {- 0, smooth wall friction factor 
FRIC } -1, no wall friction 

= roughness factor for rough walls 

BLKF = annulua effective area factor (= .83 For fully developed i 


{ 
turbine gett i 
TANS = tangent of flow separation spread angle (.1 recommended) i; 4 
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CDB = drag coefficient of struta across annulus (1-~1.2 RECM) 
SK a ratio of air specific heats, 
RG = air gas constant 
{= 0. for reverse flow annular or can combustors 
, IREAX)= 1, for axial flow annular. First data set is for OD 
| panel. Program expects a second set for ID panel 


CASE TERMINATION 
After last card of case: 


o In Column 1, Column 2 blank ~ case repeated with changes, 
next card ia title card followed by cards with changes 
from previous run. 


oo In Columns l and 2, next card is EOF to quit or new title 
card followed by all cards for complete new cage, 
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PROGRAM 117 INPUT DATA SHEET 
TOR ELEMENT VAHDS = SHEET ¢ 


INPUT FORMAT 


SEEM ALL NUMBERS MUST HAVE 
oe oe DEC TMAL POINTS 
1 


Peltor aware HOLES | NHTYe 


ee eae a ee NHTYP 


Flow passage ta divided into length 


iL) and flow (F) elements, element numbers, NELEM, 
cen be in arbitravy order, Leas, 10, Le 3. 4, 16, 30. 


The cards are stacked in order 
from inlet to last P because numbers are arbitrary, a nuw element can be inserted without 
renumbering other carda. 


bent CLEMENT 
Firet & card io annulus inlet 


aeparatuly 


for internal casee, <nd card ia dome inlet 


Le carda are used with T RISE 


SE « AT due to combuation 
in thie wlement (ne L cacda) 


for uxternal cases use only L carde 
XL, XP = X COORD to ond of element 


PLENUM 
tL = Liner 
RL, RP «= Radius to and of elument Pe Planum 
RP oo Ob, XL, RL Th) 


(Fort antermal flow Xb, 


CL © Length of element (optional) 
TLIN = Mean wall temp. over CL, °R 
tf =~ 0 then ae as 2tTi 


BLKI = Erontel Area 2) C xtruty 
Annuluse Area 
DOKI + Width of ateut 


7 y Cm 


Fo Carde ave insertud buetwoun L cards at points where flow ie oxtracted. 


ALL Dimensions in inches 


(ch = 0) = For both external and internal cases, 
Internal and external flow cases must ba run 


Li 
ery 
eee 


jp uingr atx 1 [osx + | 
[Tuer [aux 3 [ook | 


VANOTH RLIMENT 


1 @ 6% BtAaldnr 


ut) and LING 1B CALC. 


COMMON ORIGIN 
ALL SLEMENTS 


Flow conditions Lt 
into F elem. are thoae from upstream L elem, 


Types: FR « Fixed Flow Ratho,g WW) 


Fo + Fixed Orf. Diam, 
ra « 


ft « tnternal flow elem. 


| 
WeWl ie Firat Juess) 
Bleed flow (not included in liner flow) 


iy \ 
a 
(input to theee elements ie obtained from an external 
flow solution) 


WWI » Orifice Flow/Inlet Flow 
NHOLES « No. of Orificus 
NHTYP +» Hole Type (For co) 
1. Pluah Hole, Thin wall 
2. Plunged Hole 
3. Cooling skirt 
4. Fiugh Hole, Thick wal} 
$. cd Tnput 
6. Rectangular Hole 


CD @ DISCH Coefficient (NHTYP = 5 Only) 


DHOLES © Hole Diam, (for FD Only) 

I S2P = 1, Separation ie Reattached 
VIET © Jet Fiicscis, Velocity, FPS 
ANGT « Jet Injection Angle, Deg. 


ABETA » Swirl Angle in Annulue Outside Fr Elem. 


ChOPLEAAMANTR 


Wy hy 
aw 


i ba a : 


“tats 
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3D COMBUSTOR PERFORMANCE MODEL 


22 VARIABLE ‘TITLE CARDS (4A10) SAME ORDER AS IMPRINT 
1l 


CASE TITLE CARD (8A10) 


. LP], MP1 LPL IPLAX MODEL MODER IPAR ITRAD 
; 3 
IU MODEN INTAPE  IDW IRES 
4 le aad 
L SOLVE (8 (12, 8X)) 
U Vv WwW pt KE ¢ & MFU 
eee ee ee ee ee ed 
tele) ff FX FY FZ 
ICTDMA (8 (12, 8X)) 
U Vv W Pp! KE € DB MPU 
ap EODERK ERIN DORR ace EPR Gea oe RR PERRET 
MCO fF 


IMPRINT (8(IT2, 8X)) 
U Vv Ww PRESS KE £m ) MEU 


FA emi a EG isl Se et 


or ee ee ee 


eee — a eee 


‘TREMP fF Favg Ex Ly ae MCO MH20 _ 


ae ee 


MO2 MCO2 MN2 BEF DENSITY EVAP 


ELAXATION PARAMETERS (8810.4) 


U Vv Ww p! KE & 9) MEU 
eC Roe a ae [aa a CH 
MCO F FR Fy F2 PRESS DENSITY VISCOS 


LAMINAR PRANDTL NUMBERS (810.4) 


U V WwW p! KE € 9) MFU 
Ono lee ee ee ede eed 
MCO Rf 


TURBULENT PRANDTL NUMBERS (8E10.4) 


U Vv Ww p! KE € b MEU 
i a Ra Aa Dr a (Rear RT | 
MCO fF 


X-COORDINATES (1L-LP1) (8£10.4) 
ll 


ers Ea eee ors Re es eh) 
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- ewe 
~ 


‘ 
i 
| 
; 


ee eee mys S72) ae yo ee = oC fay ee! <2 Sg ab ee deed & 
d + error. 18 te C eee 
Py 
i 
+ 


l2 


Z-COORDINATES (1-NP1) (8E10,.4) 


peas aa ee ee es Ae Ra REN 


IWEL JWrIo (2 (12,8X) ) 
IWLI VALUES (8(12,8X%)) (SKIP 15 AND 16 IF IWEI = 23 


| 
| 


145 


JWLO VALUES (8 (£2, 8X)) _ 
a ON | (Ry SAE Deen Se Maa ae 


: PRESS DEN  ABSOR SCATR  AKFAC  ALFAC (BELO. 4) 
> 
: cx HY HEU FUMCO (@E10.4) 
Ie 
r. le 
e 
: PREXP1] ARCON1] CR1 PREXP2 ARCON2 CR2 BE10.4 
Mee ee 
(8E10.4) 
on C2 CD AMU ERROR ‘TCYLW TINLW TLIP 
20 
(2(13, 7X), 6(12, 8X) 
LASTEP IJUMP  JSWL  oJSW2  NUINT NVING 
21 
USW vsw SWNO  AFSW FSW TSW (BE10.4) 
22 


NFNZ LISPRAY TFUEL 


= Sea (eR (ae, RENN RIO! RS (Sar | 


XO yo 20 ALFA BETA DELTA THETA 1] THETA 2 
24 | 


RSP Wr SMD VFUEL (SKIP 24 IF NFNZ = 00) 4 


146 t 


as 


26 


27 


28 


29 


30 


32 


33 


(SKIP CARDS 25 -»- 30 IF NUINS = 00) 


I_~ LOCATION OF COOLING SLOTS aaa = 
J = LOCATION OF COOLING SLOTS = as 


AXIAL SLOT VELOCITY (8E10.4) 


ees Seles eens eee ee Eee 


TANG. SLOT VELOCITY (8E10.4) 
SLOT FLOW RATE (8E10.4) 


SLOT TEMPERATURE (8510.4) 


(SKIP CARDS 31 — 38 IF NVINI = 00) 
I - LOCATION OF RADIAL INJECTION (8 (12.8X) ) 


aca ae eee Lae rae! Ren Gee 


J - LOCATION OF RADIAL INJECTION (8(12, 8X)) 
K - LOCATION OF RADIAL INJECTION (8(12, 8X)) 


INJECTION VELOCITY (8E£10.4) 


ae ct ae ee ee 


ea, 


36 


2 ee ee ee 


INJECTION TURBULENT KINETIC ENERGY (810.4) 


Heise] Rca) eas eed ee ae ee ee 


INJECTION .TURBULENT LENGTH SCALE an eae 


INJECTION FLOW RATE (8E10.4) 


INJECTION TEMPERATURE (8E10.4) 


3D COMBUSTOR PERFORMANCE MODEL 
INPUT SHEET DESCRIPTION 


- 62 For axisymmetric geometry 4 
MODEL - 01 For laminar viscosity 


' | 

} Card Set Description 
| 1 Each card is a heading for a particular three- q 
RY} dimensional array that is printed out. These never ‘i 
in change. fa 
we a 
| 2 Case title card i 
a 3 LPI - Number of grid nodes in axial (x) dir. F 
io MPI - Number of grid nodege in radial (y) dir. }e 
i NPI - Number of grid nodes in tang. (2) dir. - a 
a IPLAX - 01 For plane geometry 4 


~ 02 For K-E viscosity model 
MODER - 01 For kinetic controlled combustion 
~ 02 For kinetic and turbulence controlled 


| 
| 
| combustion 


rm, 


IPAR - 01 For absolute pressure 
- 02 For relative pressure 
ITRAD - 01 No radiation 
~ 02 With radiation 


4 Tv ~ O1 Input units are international system Uns 
(i.e@., meters, kilograms, degrees kelvin, 
newtons, joules, radians, seconds or com- 
binations thereof) 
- 02 User selected input units 
MODEN <- 01 Density is fixed at the value of "Den" 
on Card 17 : 
~ 02 Density calculated from perfect gas law 
INTAPE ~ 00 Initial conditions not printed l 
~ 08 Initial conditions printed 
IDW - 00 Inner boundary is axis of symmetry 
- Q1 Inner boundary is wall 


IRES - 00 This is a new case 
- 01 This is a restart of previous case 


5 An 01 in proper field indicates that this particular 
variable will be solved for; an 00 indicates that 


it will not be. 


6 Indicates the number of “sweeps” made in the solve 
routine for each variable, 


eee - 
5 - - 


7 An 01 indicates that this variable will he printed, 
an 00 indicates that it will] not be, 


8 Relaxation parameters for each variable. 

9 Laminar Prandtl numbers for each variable. 
10 Turbulent Prandtl numbers for each variable. 
11 X-coordinates (LPI values) 

12 RI - Radius of inner houndary 


Y-coordinates as measured from inner 
boundary (MP1-1) values. Since Y¥(1) is 
always 00, RI is read in its place. 


| 
| 
| 
| 
| 
14 IWEI ~- I-node at which inclined wall ends 


13 Z ~ coordinates (NPI values) 
JWIO - J-node at which inclined wall starts 
15 IWLI ~ Starting I-nodea of the calculation domain 
when inclined wall is present 
16 JWLO = Ending J-nodes of the calculation domain 


when inclined wall is present 


17 PRESS ~<- System pressure 

DEN - The value of density if option MODEN = 01 
is selected 

ARSOR -=- Absorption coefficient in radiation model 

SCATR - Scattering coefficient in radiation model 

AKFAC - Internally defined turbulent kinetic 
energies are AKFAC times the appropriate 
velocity squared 

ALFAC - Internally defined turbulent length scalee 
are ALFAC times the appropriate distance 


/ 


18 Cx - Carbon atoms in fuel molecule 
| HY - Hydrogen atoms in fuel molecule 
; HFU - Heat of formation of fuel 

‘ FUMCO = Initial value assigned to Mag 
rs 19 PREXPI - Preexponent of ist reaction 


ARCONI - Activation energy divided by gas constant are 

of lst reaction (E/R) i 

CR - Constant in turbulence controlled reaction he SH 

rate for lst reaction J 4 
PREXP2 - Preexponent of 2nd reaction 
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0 ee 


20 


21 


22 


23 


24 


mepeneinaeeee, 


ARCON2 ~- Activation energy divided by gas constant 


CR2 


cl 

C2 

CB 
AMU 
ERROR 
TCYLW 
TINLW 


TLIP 


of 2nd reaction (E/R) 
- Constant in turbulence controlled reaction 
rate for 2nd reaction 


- Turbulence model constant 

- Turbulence model constant 

~ Turbulence model constant 

- The value of the viscosity if option model 
= 01 is specified. Also the laminar 
viscosity used in the "wall functions" 

~ Program will terminate if total error in 
mass becomes less than this value 

- Temperature of cylindrical portion of 
combustor wall 

- Temperature of inclined wall portion of 
combustor and of dome. 

~ Temperature of cooling slot lip. 


LASTEP - Maximum number of iterations 


IJUMP 


Jswl 
JSW2 
NUING 


NUINJ 
USW 
VSW 
SWNO 
AFSW 
FSW 
TSW 


NFN2Z 


- Number of iterations between array 
printout 

~ J-node at start of dome inlet 

- J-node at end of dome inlet 

- Number of axial injection points (cooling 
slots) 

- Number of radial injection points 


- Axial velocity of dome inlet 

- Radial velocity of dome inlet 

- Ratio of tangential to axial velocity at 
dome inlet 

- Plow rate of fuel and air through dome inlet 

- Flow rate of fuel through dome inlet 

- Temperature at dome inlet 


~ 00 No liquid fuel nozzle 
- 01 Liquid fuel nozzle present 


ISPRAY - Droplet evaporation routine is called 


TFUEL 


aod 

Yo 

2Q 
ALFA 
BETA 
DELTA 


every ISPRAY iterations 
~ Initial temperature of liquid fuel 


- X-location of origin of fuel nozzle spray 
- Y-location of origin of fuel nozzle spray 
- Z-location of origin of fuel nozzle spray 
- Nozzle cone angle 
- Nozzle back angle 
- Nozzle down angle 


THETAL - Initial spray cone segment angle 
THETA2 - Final spray cone segment angle 


tai lea ss FR aw, 


{ 
7 
\ 
RSP - Number of spray cone rays 
WE - Fuel flow rate 

SMD ~ Sauter mean diameter ‘ 
VFUEL - Initial fuel droplet velocity | 
25 I node location of cooling slots a 
26 J node location of cooling slots : 
27 Cooling slot axial velocity i 
28 Cooling slot tangential velocity ; 
29 Cooling slot mass flow rate 4 
30 Cooling slot temperature fe 
31 I node location of radial injection 
32 J node location of radial injection : 
33 K node location of radial injection } 
34 Radial injection velocity : 

35 Radial indection turbulent kinetic energy 

| 36 Radial injection turbulence length scale 
| 37 Radial injection mass flow rate 


38 Radial Injection temperature 


Te 


LINER COOLING MODEL INPUT SHERT 


CASE TITLE (12A6) 


l 
GASE TITLE (12A6) 
(Q (12, 8X) 
N KRAD MASSTR  ISWRE MODEL INERT IELUX I'l EMP 
soa cece emer SLPS) (emeSainan verona) eae ciate mae 
NS'TAT NPROF NPLOI ITEST LAS'TEP (5 (15.5X) 
3 
XU XULAS'T FRA XEND Xour PRESS POWER (8110.4) 
4 
NBP (12) 
Oly See ee ee ee ee ol ae ol 
(8h10.4) 
x RI RE x, RI RE, x, RI, 
6 Ses es 
RE x RIg RE Ke RI, RM, (8110.4) 
RA RB RD (@n10.4) 


ieee ce Neel WERE EA: evel areal See) eee ees eee 


ieee aie es (ee Be Ree eae Pee 
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iaze Uhalgs dea ok: Sawa Age oe ide Bande as mCGRaas, Gaias eat eee n Meas hee ee SN TI ae ST Da el ee | 


NAME LIST 


6 [srnpur | an as eee ae eae 


UA UD (BE10.4) 


1g CE AE maa ey EE ee Ee 


F2A F2D TA TD T (8E10.4) 
1 ced a AEE) ON ORAS NAA RGN! SNES 
NIN (T2) 


Y VALUES (GE10.4) 
l2 


U VALUES 10.4) 
Fa SS ASS TLR RTT LSE CE 
8210.4 
14 


TEMPERATURE VALUES (8E10.4) 
15 

MFU VALUES (@E10.4) 
16 . 

PHI VALUES (8E10.4) 


| 
| 
\ 
i 
H 
{ 
' 
‘ 


MCO VALUES (810.4) 


18 

TURBULENT KINETIC MNERGY (BE10.4) 
19 

TURBULENCE LENGTH SCALE (8E1.0.4) 
20 

VELA, VELAG ‘TAN, TAN oO PR, PRg (@K10.4) 


7 eee Oa KE Sane es RS 


(12, BX, 7810.4) 
Ni NZ XbpP Yop UF vi SMD wh TeUBL 


| SRS APS GANT AERO RI ap Ue ESE 


LINER COOLING MODEL 
INPUT SHEET DESCRIPTION 


Card Set Description 
1 Case title cards 


2 N - Number of cross stream intervals 

KRAD - 00 For plane geometry 

01 For axisymmetric qeometry 
MASSTR - 00 Species equations not solved 

01 Species equations solved 
ISWRL - 00 Always 
MODEL = 01 Laminar viscosity 

02 Two-equation K-E viscosity model 
INERT - 01 Species are inert 

02 Species will react 
IFLUX - 00 Wall adiabatic 

01 Wall temperature calculated 
ITEMP - 00 Isothermal 

01 Enthalpy solved 


3 NSTAT = Number of steps between printout of output 

variables 

NPROF = Number of steps between printout of output 
variables and profiles 

NPLOT = Number of steps between line-printer 
plots 

ITEST - 00 No extra printout 
01 Extra printout 

LASTEP - Maximum number of marching steps 


4 XU - Initial X-location 
XULAST - Final X-location 
FRA ~- Fraction of grid height to be used as 
step size 


XEND - X-location of end of inner wall 
XOUT - X-location of end of outer wall 
PRESS - Pressure 
POWER - Control for node spacing 
5 NBP - Number of boundary pairs | 
6 Xi - X-Location at which boundary is 
specified { 
RI, - Inner boundary radius { 
RE; - Outer boundary radius | 
- RA - Radius of axis of symmetry | 
RB ~ Inner boundary radius at initial 
X-location boa 
RD - Outer boundary radius at initial X-location i 
A | 


a i 
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ae eran 


10 


11 
12 
13 
14 
15 
16 
1? 
18 
19 


20 


21 


22 


Namelist, see page 158 for variables 


UA - 
UD - 


F2A 
F2D 
TA 

TD 
TWALL 


NIN - 
Y-values 


U-values 


Blank card(s) 

Temperature values of input initial profile 
MFU values of input initial profile 

Total fuel values of input initial profile 


M values of input initial profile 


co 


KE values of input initial profile 
(Read only iff KREAD = 1) 1 4 


fm Values of input initial profile 


SU TT TE TRE LC TL RAO eh ee 
7 i = 


Axial velocity at "free" inner boundary 
Axial velocity at "free" outer boundary 


Fuel mass fraction at "free" inner boundary 
Fuel mass fraction at "free" outer boundary 
Temperature at "free" inner boundary 
Temperature at "free" outer boundary 

Wall temperature 

Number of points on input initial profile 
of input initial profile (NIN values) 


of input initial profile 


(Read only if KREAD = 1) 4 


VELA, 
VELA, 
TAN, 
TAN, 


z 
3 
2 
Nw 
1 @f Pt 06 


< 
be} 
'eebtbpetfo 


Inner annulus velocity 
Outer annulus velocity 
Inner annulus temperature 
Outer annulus temperature 
Radius of inner plenum 
Radius of outer plenum 


00 No fuel nozzle 

01 Fuel nozzle present 
X-location of fuel nozzle 
Y-location of fuel nozzle 

Axial velocity of fuel droplets 
Radial velocity of fuel droplets 
Sauter mean diameter 

Fuel masa flow rate 

Fuel temperature 
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LINER COOLING MODEL 
NAMELIST INPUT 


; VBL Value Description 
| IUTRAP 2 Test for neaative U's, see STRIDE(2) 
ULIM 0.05 Entrainment control 
PEILIM 0.01 Max. fractional mass flow change per 
atep 
| AFAC 0.2 Relaxation on duct area deviation 
AEXDLM 0.02 Max. duct area fractional deviation per 
step 
NOVEL 2 01-U not solved for, O02-solve for U 


ARCON] 28500 Activation energy divided by gas constant 
for fuel reaction 

PREXP1 5E + 15 Preexponent for fuel reaction 

CR1 6.0 Eddy breakup constant for fuel reaction 

ARCON2 12500 Activation energy divided by gas constant 
for CO reaction 

PREXP2 6E + 8 Preexponent for CO reaction 


CR2 4.0 Eddy breakup constant for CO reaction 
cP 1048 Specific heat 
TPDDX 1 01 = Std genmix pressure grad. calculation 
02 - "Grid filling duct" version 
Cl 1.42 Turb. Constant 
c2 1.92 Turb. Constant 
CD 0.09 Turb, Constant 
C2vT 0.36 Turb. Constant 
AKFAC 0.03 Factor for internally generated kinetic 
energy profiles, KE = AKFAC*U 
ALFAC 0.02 Factor for internally generated length 
scale profiles, f/m = ALFAC*AY 
PREF -~ Turbulent Prandtl numbers 
Cx 1.0 Molecular carbon value of fuel 
HY 4.0 Molecular hydrogen value of fuel 
HFU -49317 Heat of formation of fuel ‘ 
MODER 2 1 - Kinetic only, 02 - Kinetic + Diffusion _ 
ITHIN l Thing output profiles by printing every ae 
ITHIN value ; 
| KREAN 0 0 - K & E profiles generated internally 
1 - K & fm read in a) 
ie URAT 0.83 Cooling slot velocity shape factor ty fad 
i ABSOR 0.1 Absorption coefficient in radiation model be 
EMISW 0.8 


Emissivity of liner wall 
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TRANSITION MIXING MODEL INPUT SHEET 


CASE TITLE CARD (12A6) 


1 [L 


CASE TITLE CARD oT 


(B(I2, 8X)) 


; 
| 


Sore. Pe 


ig N KRAD MASSTR ISWRL MODEL INERT IFLUX ITEMP 
4 
2 
i NSTAT  NPROF_ NPLOT _ITEST LASTEP 5 (I5, 5X 
4 3 
(8E10,4) 
ZUI ZUE ZULAST FRA  2ZEND ZOUT PRESS POWER 
4 
NBP  ICURV NRCVI NRCVE (8(12, 8X)) 
5 
(BE10.4) 
RI, XI, RI, -XIQ,° «ORI XI, RE XIg 
2 ee goal aries ROR aa ae oe eee 
(8E10, 4) 
RE, XE, — REY XE, RE, XE, REg -XEqeee 


(8E10.4) 
ay RCvI aI RCVI ai RCVI aI RCVI +6 


~s 
Aa RC mee Dae ek a Tl eS a i 


parle x. 
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wo 


10 


ll 


l2 


13 


14 


15 


16 


(8E10.4) 


RCVE l alo RCVE 4 ZE RCVE. alg ROVE 4 es 


(iene (ea a aa eer irae ieee eee 


RA RB RD (8E10.4) 


NAMELIST 


stxpur | aes ae: aera ae 


UA UD (8F10.4) 


F2A F2D Ss TA tp = Twall (8E10.4) 


cee Cae ee Peel ree eee EN 


NIN = 


Y¥ VALUES (8E10.4) 


U VALUES a a 


TEMPERATURE VALUES ae 
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a 

‘ 

Pr 
. 9 
bf 
i 


{ 
| 
ev dia shee 1th dl 


: 
| 

_ 

| gh eo) RS (SRC ORI DE Ree aa ; 
\ 4 
i 


' 

j 

4 eee, EO Eres (RN DAA Rivet aS ae - 
i 

— 

ho 

TURBULENT KINETIC ENERGY (8E1L0.4) il i 

20 


TURBULENCE LENGTH SCALE 


(8E10.4) ‘ 


: TRANSITION MIXING MODEL 
7 INPUT SHEET DESCRIPTION 


‘ Card Set Description 

K 

: 1 Case title cards 

t 2 N - Number of cross stream intervals 
Wi. 


KRAD - 00 For plane qeometry 

01 For axisymmetric geometry 
MASSTR Always 
ISWRL =~ 00 Always 
MODEL - 01 Laminar viscosity 

02 Two-equation K-F viscosity model 
INERT - 01 Always 
IFLUX - 00 Walls are adiabatic 
, 01 Wall temperature is specified 
a ITEMP - 00 Isothermal 
| 01 Enthalpy solved 


STE LEE 
t 
Qo 
Oo 


3 ’ -TAT - Number of steps between printout of output 

variables 
r | NPROF = Number of steps between printout of output 
5 variables and profiles 
NPLOT = Number of steps between line-printer 

plots 
ITEST = 00 No extra printout 

Ol Extra printout 
LASTEP - Maximum number of marching stens 


4 aut - Initial Z-location on inner houndary 
ZUE - Initial Z-location on outer boundary 


-- a ae 
PT ee eee le 


ZULAST Final Z location 
FRA Fraction of grid height to be used as 
step size 


— : =e 


ZEND - Z-location of end of inner wail 
ZOUT ~ Z-location of end of outer wall 
PRESS - Pressure 

POWER -~ Control for node spacing 


5 NBP - Number of boundary pairs 
ICURV - 00 No radius of curvature effects 
7 | 01 Radius of curvature effects 

NRCVI - Number of radius of curvature points 
j specified on inner boundary 
= NRCVE - Number of radius of curvature polnts 
= specified on outer boundary 


RI, Inner boundary radius 
XxTj Inner boundary X-location 
RE, Outer boundary radius 
XE; - Outer boundary X-location 


wee lolniee eed OS Ca ee A cee ep ee RN GR TREE Geta ee ae eat tS atm rege te tepare 
aimee Mey : 2 2 
. so ene en ore eee Leas 


7 21, - Z~location on inner boundary at which 
radius of curvature is specified 

RCVI,; - Radius of curvature of inner boundary 

zE, - Z-location on outer boundary at which 


radius of curvature is specified 


ri 

Q 

< 

a 
He 
i] 


‘ Radius of curvature of outer boundary 

‘| 

‘ 8 RA ~ Radius of axis of symmetry 

RB - Inner boundary radius at initial Z-location 
4 RD - Outer boundary radius at initial Z-location 
H 9 Namelist, see page 164 for variables 

{ 

| 10 UA ~ Axial velocity at "free" inner boundary 

‘y UD ~ Axial velocity at "free" outer boundary 


ll P2A ~ Fuel mass fraction at "free" inner boundary 
F2D ~ Fuel mass fraction at "free" outer boundary 
TA ~ Temperature at "free" inner boundary 


TD Temperature at "free" outer boundary 
TWALL Wall temperature 
12 NIN - Number of points on input initial profile 
13 Y values of input initial profile (NIN values) 
14 U values of input initial profile 
15 Blank card(s) 
16 Temperature values of input initial profile 
17 Blank card(s) 
18 Blank card(s) 
19 Blank card(s) 
4 20 KE values of input initial profile 


(Read only iff KREAD = 1) 


21 gm values of input initial profile 
(Read only if KREAD = 1) 


SOS ETRE TT 


rr CE re er 


Se, 


f 163 : 


co Ele cc it Cea ee en ee ee ee ck eer 


re 
i 
4 


VBL 


IUTRAP 
ULIM 
PEILIM 


AFAC 
AEXDLM 


NOVEL 
ARCON1 


PREXP1 
CRl 
ARCON2 


PREXP2 
CR2 
CP 
IDPDX 


erent am: -atee mee cerns = 


TRANSITION MIXING MODEL 
NAMELIST INPUT 


Description 


Test for negative U's, see STRIDE(2) 
Entrainment control 

Max. fractional mass flow change per 

step 

Relaxation on duct area deviation 

Max. duct area fractional deviation per 
step 

O01 ~ U not solved for, 02 - solve for 

U 

Activation eneryy divided by gas constant 
for fuel reaction 

Preexponent for fuel reaction 

Eddy break-up constant for fuel reaction 
Activation energy divided by gas constant 
for CO reaction 

Preexponent for CO reaction 

Eddy break-up constant for CO reaction 
Specific heat 

0). - Std genmix pressure grad. calculation 
02 = "grid filling duct” version 

Turb, Constant 

Turb. Constant 

Turb. Constant 

Turb. Constant 

Factor for internally generated turb, 
kinetic energy profiles, KE = AKFAC*U 
Factor for internally generated length 
scale profiles, fm ™ ALFAC*AY 

Turbulent Prandtl numbers 

Molecular carbon value of fuel 

Molecular hydrogen value of fuel 

Heat of combustion of fuel 

01 - Kinetic only, 02 - Kinetic + Diffusion 
Thins output profiles by printing every 
ITHIN value 

00 - KE & ?/m profiles generated internally 
Ol ~ KE & f/m read in 
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Te Pate ae 


re. eee 


CASE 'TI'TLE CARD (12A6) 


Be ee eee 


CASE TITLE CARD (12A6) 
Ee ee ee eae BE A BN 


(8(I2, 8X)) 


N KRAD MASSTR ISWRL MODEL INERT IFLUX ITEMP 
2 
NSTAT NPROF NPLOT ITEST DLASTEP fs 5 (15 5X 
Se SE OR] a TS 
XU XULAST FRA XENDC xourT PRESS POWER BE10.4 
4 


RA RB RD 8E10.4 


¥ NAMELIST 


j a [srupur| | ren Bee Be ae 


- | UA UD VTA _VTD 8E10.4) 


: F2A F2D TA TD  ‘Twall (BE10.4) Ll 

_ : 
NIN NUI NUE NVI NVE (5(12, 8X)) . 

11 


Y VALUES (8E10.4) 
12 


U VALUES (BE1LO.4) 
13 


V§ VALUES (8E10.4) 


14 


x TEMPERATURE VALUES (810.4) 
15 


g MEU_VALUES (8E10.4) 


a 16 
MCO, VALUES (BE10.4) 


| 


q 166 | 
{ | 


o ate te Seeen aernmemepmgetgeners terete = 


MCO VALUES (8E10.4) 


{ 18 
MOX VALUES (8ELO.4) 
19 
MH,O VALUES (BE10.4) 
1 Rae (ares GE CRO Dann ane IMRT Rie 
MH, VALUES (8B10.4) 
21 esis a eR ae I Gene a ES 
TURBULENT KINETIC ENERGY (BE10.4) 
22 
TURBULENCE LENGTH SCALE (SELO.4) 
23 


SKIP FOLLOWING CARD SET IF NUI = 0  cemmewmsmmcmmnenenmee 
X - LOC, OF INTERNAL COOLING SLOTS (8E10.4) 


| a SAS EN RUS Ieee re elie See Re 


LIP LENGTH OF INTERNAL COOLING SLOTS (8E10.4) 
25 
U_- VELOCITY OF INTERNAL COOLING SLOTS (8E10.4) 
26 
VT - VELOCITY OF INTERNAL COOLING SLOTS (8E10.4) 
27 
167 


| 


me oe 
pei Eeoeh mete 


TEMPERATURE OF INTERNAL COOLING SLOTS (BE10.4) 


FLOW RATE OF INTERNAL COOLING SLOTS (8E10.4) 
29 

SLOT HEIGHT OF INTERNAL COOLING SLOTS (8E10.4) 
30 

SLOT TO METERING AREA RATIO FOR INT SLOTS (8E10.4) 
31 


SKIP FOLLOWING CARD SET IF NUE = Q  aeenmmevmemmcnecs 


X-LOC_OF EXTERNAL COOLING SLOTS (BE10. 4) 
32 _ 

LIP LENGTH OF EXTERNAL COOLING SLOTS __(BE10.4) 
33 ; 

U - VELOCITY OF EXTERNAL COOLING SLOWS (BEL0.4) 
34 

V,, ~ VELOCITY OF EXTERNAL COOLING SLOTS (BE10.4) 


St ee 


(aes Sia ea RR ae ae Rene ieee Dee) 


t 


ee 
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ae Ee PER) 


' 
t 
5 een wire ‘ oY 
tec cet eae Rane ne Mn ctl i btn dade anne ta B tena le a ee a tle i el 


Re PETER op reve voce 


i, oe le epee . = 
. oe aaa ~~ 
‘ es: 16) 


. TEMPERATURE OF EXTERNAL COOLING SLOTS 8E10.4 
36| 

FLOW RATE OF EXTERNAL COOLING SLOTS sas uae 
37 


SLOT HEIGHT OF EXTERNAL COOLING SLOTS (8E10,4) 
38 

SLOT TO METERING AREA RATIO FOR EXT SLOTS (8E1L0. 4) 
39 

SKIP FOLLOWING CARD SE. IF NVI =~ 0 

X ~ LOC OF INTERNAL RADIAL INJECTION (8E10.4) 
40 

U = VELOCITY OF INTERNAL RADIAL INJECTION (8E10.4) 

V_- VELOCITY OF INTERNAL RADIAL INJECTION (BE10.4) 
42 

VT = VELOCITY OF INTERNAL RADIAL INJECTION (8E10.4) 


Oe ee 
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TEMPERATURE OF INTERNAL RADIAL INJECTION (8E10.4) 


FLOW RATE OF INTERNAL RADIAL INJECTION (BE10.4) 
45 


SKIP FOLLOWING CARD SET IF NVE = QO 
X_=- LOC. OF EXTERNAL RADIAL INJECTION (BE1L0.4) 


js 
ate Be eh Stan anabarne oa 


46 i 
3 
U = VELOCITY OF EXTERNAL RADIAL INJECTION (8E10.4) | ; 
47 
V = VELOCITY OF EXTERNAL RADIA INJECTION (BELO.4) 1 
48 
Li 
VT - VELOCITY OF EXTERNAL RADIAL INJECTION (BE10.4) | 
49 | 
| 
TEMPERATURE OF EXT'ERNAL RADIAL INJECTION BE10.4) | 
50 
FLOW RATE OF EXTERNAL RADIAL INJECTION (8E10.4) 


51 


EMISSIONS MODEL 
INPUT SHEET DESCRIPTION 


{ Card Set Description 
1 Case title cards 
2 N - Number of cross stream intervals 
KRAD - 00 For plane geometry 
i 01 For axisymmetric geometry 


MASSTR 


t 
ro} 
jo] 


Species equations not solved 

01 Species equations solved 
ISWRL - 00 Swirl velocity not solved 

01 Swirl velocity solved 
MODEL - 01 Laminar viscosity 

02 Two-equation, K-E, viscosity model 
INERT - 01 Species are inert 

02 Species will react 
IFLUX <- 00 Wall adiabatic 

01 Wall temperature specified 
ITEMP = 00 Isothermal 

01 Enthalpy solved 


3 NSTAT - Number of steps between printout of output 
variables 
| NPROF - Number of steps between printout of output 
variables and profiles 
| NPLOT ~ Number of steps between line printer 
plots | 
| ITEST = 00 No extra printout | 
O1 Extra printout 
| LASTEP ~- Maximum number of marching steps 
ay | 
| 


4 XU - Initial X location 
XULAST - Final X location ; 
PRA - Fraction of grid height to be used as : 
step size 
XEND X location of end of inner wall 


{ | XxXOUT - X¥ location of end of outer wall 
§ PRESS = Pressure 


POWER Control for node spacing 
5 NBP - Number of boundary pairs 
6 X - X location at which boundary is specified 
Ri, - Inner boundary radius 
RE; - Outer boundary radius 
7 RA ~ Radius of axis of symmetry 
, RB - Inner boundary radius at initial X-location 


RD - Outer boundary radius at initial X-location 


8 Namelist, see page 174 for variables 


9 UA - Axial velocity at "free" inner boundary 
UD - Axial velocity at "free" outer boundary 
VTA - Tang. velocity at "free" inner boundary 
VTD - Tang. velocity at "free" outer boundary 

10 F2A - Fuel mass fraction at “free” inner boundary 
F2D - Fuel mass fraction at "free" outer boundary 
TA - Temperature at “free" inner boundary 
TD - Temperature at "free" outer boundary 
TWALL =- Wall temperature 

11 NIN - Number of points on the input initial 

profile 

NUT - Number of cooling slots on inner boundary 
NUE = Number of cooling slots on outer boundary 


NVI - Number of radial injection points on 
inner boundary 

NVE ~ Number of radial injection points on 
outer boundary 


12 Y values of input initial profile (NIN values) 
13 U values of input initial profile 

14 Vg values of input initial profile 

15 Temperature values of input initial profile 
16 Mey values of input initial profile 

17 Mao2 values of input initial profile 

18 Mao values of input initial profile 

19 Mox values of input initial profile 

20 My20 values of input initial profile 

21 My values of input initial profile 

22 KE values of input initial profile 


(Read only if KREAD = 1) 


23 ém values of input initial profile 
(Read only iff KREAD * 1) 


a 24 to 31 Information pertaining to cooling slots on inner 
boundary 
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én wie 
t. 
t 


32 to 39 Information pertaining to cooling slots. on outer 


boundary 
40 to 45 Information pertaining to radial injections on 
inner boundary 
46 to 51 Information pertaining to radial injections on 
outer boundary : 
i 
! 
; i 
| | 
l 
f 
1 
! ‘ 
(3 
i 
: 
| 
4 y 
ni 
i" 
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e A 
. if 
r e Os: ar: To On tue BUT ry oh OO 


Pe eee PAR 


! EMISSIONS MODEL 
NAMELISY INPUT 


| VBL Value Deseription 
\ IUTRAP 2 Test for negative U's, see STRIDE(2) 
: ULIM 0.05 Entrainment control 
PEILIM 0.01 Max. fractional mass flow change per 
step 
: AFAC 0.2 Relaxation on duct area deviation 
: AEXDLM 0.02 Max. duct area fractional deviation per 
step 
NOVEL 2 O01 - U not solved for, 02 = Solve for 


U 
ARCON1 28500 Activation energy for fuel reaction 
PREXP1 5E + 15 Preexponent for fuel reaction 


CR1 6.0 Eddy breakup constant for fuel reaction 

CP 1048 Specific heat 

IPDDX 1 01 = Std genmix pressure grad. calculation 
02 - "grid filling duct" version 

Cl 1,42 Turb. Constant 

C2 1.92 Turb. Constant 

CD 0.09 Turb. Constant 

CovtT 0.36 Turb. Constant 

AKFAC 0.03 Factor for internally generated kinetic 
energy profiles, KE = AKFAC*U 

ALFAC 0.02 Factor for internally generated length 
scale profiles, é/m = ALFAC*aAy 

PREF -- Turbulent Prandtl numbers 

Cx 1.0 Molecular carbon value of fuel 

HY 4.0 Molecular hydrogen value of fuel 

HFU ~49317 Heat of formation of fuel 

MODER 2 l =~ Kinetic only, 02 - Kinetic + Diffusion 

ITHIN 1 Thins output profiles by printing every 
ITHIN value 

KREAD 0 0 - K & E profiles generated internally 
1- K & fm read in 

URAT 0.83 Cooling slot velocity shape factor 

TERM) 0.1 Control on specie equation step size 

TERM 2 1.EB=-4 Control on specie equation satep size 

ISMAX 500 Maximum specie equation steps 

EFU 1.0 Power on fuel in fuel reaction rate 

Fo 0.5 Power on Oo in fuel reaction rate 

Ey 0 0.5 Power on H30 in fuel reaction rate 

2 
ERO 2.0 Power on density in fuel reaction rate 
EPR 0.0 Power on pressure in fuel reaction rate 
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FUEL INSERTION MODEL INPUT 


j TITLE 80 


AIR 
PRIM SEC. ASSIST 
i we FLOW FLOW CONE SHROUD PRIM SEC 


ATOM AIR NO. NO. ANGLE EFF ORIFICE ORIFICE 
TYPE ASSIST(JP4) (JP4) DEG AREA DIA, DIA, 


Bru/ el IN? IN IN 


l~-6—~ 1) 31 41 51 61 5 80 
2 


*ATOMIZER TYPE: 00001 = SIMPLEX 
00002 = DUAL ORF 
00003 = AIR BLAST 
*®AIR ASSIST: 00001 = NO ASSIST 
00002 * WITH ASSIST 


*FOR ATOM TYPE = 00002 (DUAL ORIF) ONLY (LEAVE OUT FOR OTHERS) 
INPUT SECONDARY FLOW SCHEDULE 


Ws = SECONDARY FUEL FLOW, LB/HR 


APs ™ APsec, orntr * “Ppyow piviper vave! P/D 
CRACK 
POINT 
FLOW 
w W W W W. 
82 
yh. 44 oy. 85a OF aa OP ay 
of Gee OREN] ON Ne SNE) EME aL 
CRACK 
PRESSURE eae Uae 
Rs: iaaeer) ee eee? Mee | Sea 


*FOR ATOM TYPE = 00003 (AIRBLAST) ONLY (LEAVE OUT FOR OTHERS) 
FILMING NOZZLE FLOW AIR 


DIA ey) RFLOW Ow TEMP 
IN. LB/S me IN °R 
i 1i 31 
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Ldiadiatial AIR AIR 


*** AIR FUEL FUEL FUEL ASSIST ASSIST 
FUEL FLOW TEMP FLOW AP APs TEMP 
TYPEOPTION °R  LB/HR PSI PSI °R 
1 6 11 21 31 41 51 61 7 
6 
Rw eeKE 
ne tes ed ae 0.00001=UNIFORM GAS STREAM 
cprton| 0+0000282-D FIELD OPTION 
t-] ry 
Tongs Gas Pansy.) MAK OMAK, a 
1 11 FPS 21PSIA 31 IN. 4 °R Sl 
7 Del 
i i i 
%yoz ¢ ‘wo ¢ Pans “wax ‘wax “So 2-2 FIELD OPTION 
IN. IN.  PSIA IN, IN. 
CARDS 8 THROUGH 13 SKIPPED IF AIR FLOW OPTION = 00001 
IN = NO, OF X-DIR POINTS IN 2-D FIELI 
IN IN 2(12, 8X JN = NO. OF Y= DIR POINTS 
B 


X VALUES (8E10.4) X - LOCATIONS OF GRID POINTS (FT) 


2 Reco ici ECan URS RRS ISR Doe 


Y VALUES (B8E10.4) Y - LOCATIONS OF GRID POINTS (FT) 
10 


READ ALONG + X LINES STARTING WITH i 


U VALUES (BE10.4) SMALLEST Y VALUE (FT/SEC) 
11 | (een anata Renal Sa teal oer ea 
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Mt BC ae RT ae et a a a Ser ee re a 


; 


READ ALONG + X LINES STARTING WITH 
VY VALUES (8E10.4) SMALLEST Y VALUE (FT/SEC) 


Maen) Garis rei ence eae eee aes See 


eed REA ee Meee Sees ae Eo ee 


READ ALONG + X LINES STARTING WITH 7 
T VALUES (8E10.4) SMALLEST Y VALUE (DEG. R) ae 


13 


abd ote a aati sh NT a ata ir 


ed ta La eaten 
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APPENDIX B 


LISTING OF ANNULUS LOSS MODEL 
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